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FOREWORD 


The  remote  sensing  technique  which  is  discussed  in  this  report  represents  a  new  capability  which  could 
potentially  have  several  important  applications.  The  use  of  the  term  "remote  sensing"  we  believe  is  appropriate 
since  a  sensor  unit,  located  at  a  known  depth  in  the  ocean  is  used  (during  daylight  hours)  to  determine  the  optical 
attenuation  properties  of  the  water  above  the  sensor  and  of  the  atmosphere  above  the  water.  Knowing  the  time,  date, 
latitude,  longitude,  and  the  extra-terrestrial  solar  spectral  irradiance,  (all  easily  obtained)  it  is  a  straightforward  matter 
to  compute  the  total  path  transmittance  at  a  single  wavelength,  i.e.  a  single  value  representing  the  combined  water 
and  atmospheric  attenuation.  The  problem  of  determining  the  two  component  parts  of  the  attenuation  separately  was 
the  challenge. 

Our  approach  was  to  make  measurements  at  two  or  more  wavelengths  and  use  these  measurements  in  a  set  of 
equations  which  can  be  solved  for  the  separate  attenuations.  This  required  models  describing  how  the  attenuation 
coefficients  of  sea  water  and  of  the  atmosphere  varied  with  wavelength  for  various  turbidities.  Atmospheric  models 
existed  which  had  been  used  successfully  by  the  Visibility  Laboratory  and  others  to  infer  and  remove  the  effects  of 
the  atmosphere  on  the  water-leaving  radiances  arriving  at  the  Coastal  Zone  Color  Scanner  satelite  sensor.  Simplified 
versions  of  these  atmospheric  models  worked  satisfactorily.  The  problem  was  that  there  was  no  existing  comparable 
analytic  model  for  seawater.  The  often  quoted  and  used  water  types  of  Jerlov  are  described  by  a  table  of  attenuation 
values  specified  at  2Snm  intervals  throughout  the  visible  spectrum  for  S  oceanic  and  S  coastal  water  types.  To 
devise  a  model  around  the  Jerlov  table  would  have  required  fitting  analytic  functions  to  his  attenuation  values, 
putting  the  model  another  step  removed  from  the  actual  data  set  upon  which  the  table  was  based.  Furthermore  we 
had  available  to  us  a  large  and  high  quality  data  base  of  spectral  irradiance  profiles  which  had  been  obtained  by 
investigators  at  the  Visibility  Laboratory,  in  France  at  Prof.  Morel's  laboratory,  and  by  Okami  and  his  colleagues  in 
Japan.  These  data,  which  were  obtained  in  a  wide  variety  of  oceanic  and  coastal  waters,  were  acquired  at  wavelength 
intervals  of  approximately  5  nm  thus  permitting  a  better  representation  of  the  shape  of  the  spectral  attenuation 
coefficient  in  the  resulting  model.  That  model,  described  in  Austin  and  Petzold,  1984  (Ref.  2  and  Appendix  2, 
attached),  allows  the  complete  specification  of  the  diffuse  attenuation  coefficient  at  all  wavelengths  in  the  visible 
spectrum,  given  the  value  of  the  coefficient  at  any  one  wavelength  in  that  spectral  region. 

The  technique  which  was  developed  uses  that  analytic  model  for  the  spectral  attenuation  of  the  water  and  the 
results  have  been  very  promising  indeed.  This  report  describes  the  technique  and  examines  the  errors  that  might  be 
expected  due  to  errors  in  radiometric  calibration,  wavelength  assignment,  depth,  solar  zenith  angle,  and  in  solar  and 
atmospheric  constants.  It  also  addresses  the  matter  of  optimum  wavelength  selection.  A  subsequent  report  will 
provide  an  evaluation  of  the  technique  using  data  sets  obtained  from  various  cruises  on  surface  ships.  Among  the 
applications  we  see  for  this  remote  sensing  technique  are  the  following. 

•  The  determination  of  the  amount  of  attenuation  above  a  submerged  sensor  from  a  set  of  measurements 
made  at  a  single  point  in  time  and  space.  Separate  determinations  of  the  amount  of  attenuation  due  to  the 
water  and  due  to  the  atmosphere  are  provided. 

•  A  rather  precise  assessment  of  the  mean  diffuse  attenuation  coefficient,  K,  from  the  ocean  surface  to  the 
depth  of  the  sensor  can  be  obtained  on  a  continuous  basis. 

•  The  technique  is  particularly  useful  for  applications  involving  installations  on  towed  sensor  systems  and 
on  unattended  buoys  or  bottom  mounted  sensor  systems. 

For  the  genesis  of  this  technique  we  must  credit  Dr.  Matthew  B.  White  of  the  Office  of  Naval  Research.  It  was 
he  who  asked  the  critical  question:  "Isn't  there  some  way  you  can  determine,  using  a  submerged  sensor,  the  separate 
losses  due  to  the  water  and  to  the  atmosphere  without  moving  the  sensor  vertically?  Perhaps  use  more  than  one 
wavelength."  or  words  to  that  effect  With  that  challenge  we  have  proceeded  to  develop  the  Submerged  Remote 
Sensing  (SRS)  technique  together  with  the  required  water  attenuation  model,  both  of  which  potentially  have 
important  future  applications. 
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ABSTRACT 

Prior  analysis  and  experimentation  has  provided  strong  support  for  determining  the  attenuation  of  optical  radiation  for  the 
atmosphere  and  the  water  column  above  a  submerged  platform  by  measuring  the  absolute  down  welling  irradiance  at  two 
wavelengths.  The  technique  requires  knowing  the  two  irradiance  values  and  in  addition  (1)  the  spectral  irradiance  of  the  sun 
outside  the  atmosphere,  (2)  the  solar  zenith  angle  and  (3)  the  depth  at  which  the  irradiances  are  measured. 

With  this  capability  in  hand  the  questions  considered  in  this  report  are  which  wavelengths  to  use  and  the  effect  of  errors  in  the 
measurement  of  irradiance  and  other  parameters. 

1.0  INTRODUCTION 

The  concept  and  theory  for  the  determination  of  the  optical  properties  of  the  atmosphere  above  the  ocean's  surface  from 
spectral  measurements  of  the  natural  light  made  below  the  surface  was  presented  in  detail  in  Ref.  1*. 

One  requirement  to  achieve  a  solution  is  that  there  be  a  known  spectral  relationship  between  the  diffuse  attenuation 
coefficient,  K(X)  for  at  least  two  wavelengths  which  is  valid  for  the  various  "types''  of  natural  waters  found  in  open  oceans. 
Preferably,  a  way  of  relating  K(X)  for  all  wavelengths  over  the  range  of  interest  would  allow  flexibility  in  the  choice  of 
wavelengths  at  which  to  make  the  measurement  From  the  results,  the  K(X)  at  any  other  wavelength,  within  the  range 
covered,  could  be  determined.  A  study  of  a  large  body  of  spectral  K(X)  data  obtained  by  several  investigators  has  yielded  a 
model  which  can  be  used  for  these  purposes.  The  model  and  it's  basis  are  presented  in  Ref.  2. 

The  scheduling  of  a  field  experiment  prompted  asking  the  question  of  what  wavelengths  to  select  Five  photometers  were 
available  for  this  experiment  and  as  a  result  five  wavelengths  were  selected.  Unfortunately,  there  is  not  a  straightforward 
answer  that  is  "best"  for  all  situations.  The  "type"  of  water  is  critical  to  the  selection  of  wavelengths  and  since  this  is 
unknown  the  selection  must  be  based  on  predictions  of  the  type  of  water  expected.  The  effect  of  errors,  particularly 
photometric,  has  played  a  large  role  in  the  wavelength  selection.  If  the  photometric  measurements  could  be  made  without 
error,  then  the  criteria  for  wavelength  selections  could  be,  by  and  large,  those  that  would  allow  useable  results  to  be  acquired 
at  maximum  depth.  Perversely,  the  use  of  wavelengths  which  would  allow  the  system  to  perform  at  greater  depths  tend  to 
make  the  system  very  sensitive  to  photometric  error.  As  we  gain  more  field  experience  and  develop  more  precise  and  stable 
radiometers  for  this  application,  we  can  expect  to  obtain  greater  depth  capability. 
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2.0  METHOD 


Using  our  knowledge  of  the  optical  characteristics  of  the  atmosphere  and  seawater  it  is  possible  to  compute,  for  a  given 

wavelength,  X,  the  transmission  of  light  downward  throi*"h  the  atmosphere,  the  air-water  interface,  and  through  the  water  to 
some  depth,  z,  for  a  given  set  of  conditions.  Using  known  values  for  the  spectral  irradiance  from  the  sun  outside  the 
atmosphere,  E0(X),  we  can  obtain  an  estimated  value  for  the  downwelling  irradiance  at  depth  z,  EZ(X),  for  any  wavelength  for 
the  given  set  of  conditions.  This  is  done  at  two  wavelengths  to  produce  a  pair  of  irradiance  values  Ez(Xi),  Ej(Xi).  Then,  using 
only  these  values,  the  procedure  given  in  Refs.  1  &  2  is  used  to  perform  the  inverse  computation  to  calculate  the  diffuse 
attenuation  of  the  water,  K(X),  and  the  optical  depth  of  (or  transmission  through)  the  atmosphere  at  either  of  the  wavelengths. 
These  computed  attenuation  values  are  then  compared  with  the  original  given  or  "true  values"  to  get  an  estimate  of  the 
magnitude  of  error.  This  is  the  error  which  exists  for  a  perfect  system  and  results  from  assuming  for  the  inverse  computation 
that  the  optical  properties  of  the  aerosols  in  the  atmosphere  are  not  wavelength  dependent  This  assumption  is  necessary  to 
enable  the  solution  from  measurements  at  two  wavelengths. 

Other  sources  of  error  will  be  present  in  a  real  system  due  to  such  causes  as  radiometric  accuracy,  errors  in  the  values  used  for 
the  solar  input,  solar  zenith  angle,  depth  measurement,  etc.  Offsets  commensurate  with  those  types  of  errors  have  been 
inserted  into  the  calculations  to  find  the  sensitivity  of  the  system  to  each  of  those  sources  of  error. 

Note:  In  the  figures  which  follow  the  "true"  values  are  indicated  by  a  straight  solid  line,  the  computed  values  for  a  perfect 
two  wavelength  system  are  shown  with  a  diamond  (0)  symbol,  and  the  results  of  applying  an  error  to  some  particular 
parameter  are  plotted  using  (+). 


2.1  The  Forward  Calculation 

The  two  values  for  E^XO  and  are  determined  as  follows: 

For  values  of  t,  (X)  <1  the  transmittance  for  irradiance  through  the  atmosphere  is  taken  to  be 


Ta  (X)  »  exp  -  [.48t„(X)  ♦  To(X)  +  0.5(l-g)ta(X)]/Mo 


M) 


where 


X  indicates  wavelength 


Ta(X)  is  the  total  path  transmission  through  the  atmosphere 
Mo  is  the  cosine  of  the  sun's  zenith  angle 


Tr(X),  t„(X),  ta  (X)  are  the  optical  depths  for  unit  air  mass  due  to: 

air  molecules  (Rayleigh  scattering), 
ozone  absorption,  ana  aerosols  respectively 


g  is  an  asymmetry  factor  for  aerosol  particles  and  is 
approximately  equal  to  2/3. 


The  transmittance  for  irradiance  through  the  water  column  is 


where 


TW(X)  *  e-K<x>* 


K(X)  is  the  diffuse  attenuation,  and 
z  is  depth  of  water  over  sensor. 


(2) 


A  value  of  0.98  is  used  for  the  irradinace  transmittance  through  the  air-water  interface  for  all  wavelengths.  While  this  is  not 
strictly  true,  any  variation  from  this  value  will  not  be  large  and  would  have  a  very  small  effect  upon  the  final  results. 

The  inadiance  at  depth  z  is  then: 


E*(X)  -  0.98  •  TA(X)  ■  TW(X)  •  Eq(X)  . 


(3) 


where 


E0(X)  is  the  solar  irradiance  outside 

the  atmosphere  at  a  nominal  wavelength  X 
and  over  the  band  width  of  the  sensor. 

In  this  simulation  two  wavelengths  are  selected  and  a  spectral  bandwidth  is  chosen  for  the  radiometric  measurements.  The 
radiometric  response  is  taken  to  be  uniform  within  this  band  width.  In  a  real  system  the  spectral  response  of  the  radiometers 
will  have  to  be  carefully  determined  and  the  mean,  or  effective  wavelengths  calculated.  Also  it  will  be  necessary  to  calculate 
the  proper  effective  values  for  the  spectrally  dependent  parameters  used  such  as  EjX). 

The  value  for  EjX)  was  obtained  from  the  work  of  Neckel  and  Labs  Ref.  4.  Figure  1  is  a  plot  of  these  data.  The  data 
fluctuate  markedly  with  wavelength  and  can  not  be  represented  by  simple  curve  fitting  procedures.  Therefore,  E0(X)  is 
calculated  by  interpolation  at  one  nanometer  increments  over  the  passband  and  the  average  value  used. 

The  Rayleigh  (molecular)  optical  depth,  tr(X),  for  diffuse  propagation  of  light  through  the  atmosphere  varies  very  nearly  with 
the  inverse  fourth  power  of  wavelength  and  is  calculated  using 


x„(X) .  0.0A40J610)-4  . 


The  effective  value,  xR(X),  over  the  pass  band  Xt  to  Xj  is: 


tR(X)-0.044(Xl-Xl)1 


J  (X/670)4  •  dX  -  2.956  •  10*  (X^-X,)'1  •  (X?  -  ) 

h 


(4) 


The  work  of  Klenk  et  al.  Ref.  3.  is  used  to  find  the  spectral  absorption  coefficicent  for  ozone  as  a  function  of  latitude  and 
season.  Figure  2  shows  the  spectral  data  (+)  and  a  function  curve  which  adequately  describes  the  absorption  per  unit 
atmosphere  -  centimeter,  (  a/aun-cm).  The  total  spectral  absorption  coefficient  may  be  found  by  multiplying  these  values  by 
"total  ozone"  which  is  a  function  of  latitude  and  season.  A  value  of  0.300  atm-cm  was  used  as  representative  for  mid¬ 
latitudes.  The  optical  depth  t0(X)  is  taken  to  be  equal  to  the  total  absorption  coefficient  and  may  be  calculated  from: 


Xq(X)  -  0.300  •  0.13879  •  e*AW 

(5) 

A  -  0.0014717  ■  (|X  -  589.75I)1-5301. 

(6) 

Where  the  effective  value,  T,  over  the  pass  band  Xj  to  Xj  is 

Tq(X)  -  (X2  -  Xj)"1  •  J  'Cq(X)  •  dX 


(7) 
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The  summation  is  done  at  one  nanometer  increments  over  the  bandwidth  and  the  average  value  used  for  x0(X)  . 

The  optical  depth  for  the  aerosol,  xa(X)  has  a  spectral  dependency  which  appears  to  follow  the  power  law 

x.(X,)  =  (X./X^a  •  .  (8) 

The  exponent,  a,  known  as  the  Angstrom  exponent,  is  dependent  upon  the  amount  and  quality  of  the  aerosol.  It  falls  in  the 
range  of  1.1  to  1.4  for  clear  atmospheres  and,  approaches  zero  for  very  large  aerosols  (fog,  overcast). 


Two  types  of  atmosphere  are  used  for  illustration  purposes,  as  follows: 


Number 

Designation 

Type 

Designation 

ta(490) 

a 

1 

2 

Haze 

Heavy  Overcast 

0.10 

10-20 

1.0 

0 

For  the  "Heavy  Overcast”  case,  the  term  0.5  (1-g)  x,(X)  in  Eq.  (1)  for  transmittance  (see  Section  2.1)  does  not  hold  true.  The 
transmittance  of  diffuse  light  is  nearly  spectrally  independant  for  a  heavy  overcast  consisting  of  fog  or  complete  cloud  cover. 
For  this  case  the  transmittance  of  diffuse  light  through  the  aerosol  is  taken  to  be  10%  and  the  constant  2.3  i.e.  -In  0.1,  is 
substituted  for  the  term  involving  x,(X). 

In  some  of  the  figures  the  number(s)  appearing  on  the  right  side  are  the  number  designation  for  the  type  of  aerosol  for  which 
the  information  displayed  applies. 

The  calculations  used  to  obtain  the  irradiance  values  at  the  two  wavelengths  at  depth  z,  viz.  E^X])  and  E^fXj)  and  then  using 
only  these  two  values  to  estimate  the  atmosphere  and  water  properties,  require  that  a  relationship  exist  between  the  diffuse 
attenuation  coefficients,  K(X),  of  the  water  at  these  two  wavelengths.  For  the  forward  calculation,  i.e.,  propagating  the  light 
downward,  a  "type"  of  water  must  be  selected.  The  value  of  K(490)  is  used  to  specify  the  water  type  and  from  this,  values  for 
K  (X,)  and  K  (X^)  may  be  calculated  using  the  model  and  procedure  presented  in  Ref.  2. 

2.2  The  Inverse  Calculation 

The  model  from  Ref.  2  is  also  used  to  obtain  the  function  relating  K(Xj)  and  K(X2)  required  for  the  inverse  solution  to  the 
problem,  given  only  E^X^  and  EjfXj). 

Neither  of  the  wavelengths  X!  or  X2  at  which  the  two  irradiance  measurements  are  made  needs  to  be  at  the  wavelength  for 
which  we  would  like  to  know  the  atmospheric  and  water  optical  properties.  To  estimate  the  diffuse  transmission  of  light 
through  the  atmosphere  at  some  other  wavelength,  X3,  we  may  use  the  same  approximate  assumption  used  to  enable  the 
inverse  solution  of  the  problem  which  is  that  the  optical  depth  of  the  aerosol,  x,(X),  is  independent  of  wavelength,  hence, 
t,(X,)  -  x.fXj)  -  x,(X]).  In  effect  the  calculated  xa(X)  applies  equally  well  at  any  wavelength.  The  error,  caused  by  this 
assumption  is  not  large  enough  to  preclude  the  usefulness  of  the  procedure.  Using  Eq.  (4)  and  (5)  we  can  calculate  xR(X3)  and 
x„(X]).  From  Eq.  (1)  for  diffuse  atmospheric  transmittance  and  using  the  approximation  that  the  aerosol  term  is  constant  we 
can  derive: 

Ta(Xj)  -  Ta(X,)  •  exp  -  {  .48[xr(Xj)  -  xR(X,)]  +  x^Xj)  -  x„(X)}  (9) 

Thus  we  can  obtain  an  estimate  for  the  vertical  diffuse  transmittance  through  the  total  atmosphere  at  wavelengths  other  than 
the  measurement  wavelengths.  For  the  water  we  can  find  K(X3)  from  K(Xt)  or  KfXj)  using  the  model  for  K(X). 

3.0  CONSIDERATIONS  IN  THE  SELECTION  OF  WAVELENGTHS 

Two  criteria  were  considered  in  making  the  wavelength  selection.  The  first  and  most  obvious  was  to  achieve  maximum  depth 
at  which  the  radiometric  measurements  can  be  made.  The  second  was  to  minimize  the  sensitivity  to  system  error.  The  best 


pair  of  wavelengths  to  use  is  highly  dependent  on  the  type  of  water  for  both  of  these  criteria.  The  wavelengths  which  might 
be  chosen  for  use  in  clear  waters  would  not  be  a  good  choice  for  more  turbid  waters.  Also  these  two  criteria  are  not 
compatible.  Using  wavelengths  to  achieve  maximum  depth  capability  can  result  in  the  system  being  very  sensitive  to 
radiometic  error,  particularly  those  induced  by  wavelength  calibration,  and  using  wavelengths  which  minimize  sensitivity  to 
system  error  will  lessen  the  depth  at  which  the  radiometric  measurement  can  be  made. 

It  is  felt,  at  this  time,  that  the  second  criteria  should  be  dominant  in  the  wavelength  selection,  i.e,  do  the  best  we  can  to 
enhance  the  probability  of  acheiving  good  results.  After  it  has  been  demonstrated  that  the  radiometric  measurements  can  be 
made  with  the  precision  required  to  obtain  satisfactory  results  over  the  long  path  lengths  and  high  attenuations  involved,  then 
we  can  work  toward  greater  depth  capability. 

To  illustrate  the  effect  of  wavelength  selection  on  sensitivity  to  radiometric  error,  water  with  a  diffuse  attenuation  coefficient 
of  K(490)  =*  0.067  (Jerlov  Type  II)  will  be  used.  If  we  are  trying  for  greatest  depth  capability  we  should  choose  wavelengths 
in  the  region  where  light  penetration  is  the  greatest.  Figure  3  shows  typical  values  of  K(X)  plotted  as  a  function  of 
wavelength,  X ,  for  seven  Jerlov  water  types.  Figure  4  shows  the  "optimum"  wavelength,  i.e.,  the  wavelength  at  which  the 
minimum  K(X)  occurs,  for  water  types  in  the  region  from  pure  sea  water,  K(490)  =  0.022,  to  K(490)  »  0.250.  Considering 
only  the  water  transmission  properties  we  find  that  for  water  with  a  K(490)  «  0.067  the  "optimum"  wavelength  is  in  the 
region  of  489nm  and  we  would  select  two  wavelengths  on  either  side  of  this  value.  Figure  S  gives,  as  a  function  of 
wavelength,  the  depth  at  which  the  downwelling  irradiance  is  0.0002|iW  ■  cm*2  -  nm*1*. 

This  somewhat  arbitrary  number  allows  depth  limit  computations  which  are  useful  in  wavelength  selection.  Ideally,  in  a 
"balanced”  system  both  photometers  would  run  out  of  adequate  signal  at  the  same  depth.  Figure  5  helps  in  the  selection  of 
two  wavelengths  in  an  attempt  to  acheive  this,  for  water  with  a  K(490)  -  0.067.  The  solar  input  outside  the  atmosphere, 
E0(X),  and  the  diffuse  transmittance  through  the  atmosphere  are  included  in  the  computations  used  for  Fig.  S.  The  spectral 
solar  irradiance,  E0(X),  and  the  effective  spectral  diffuse  optical  depths,  tR(X)  (Rayleigh),  T0(ozone),  and  T,(aerosol)  are  all,  in 
this  case,  computed  over  a  lOnm  pass  band.  Also  a  clear  air  case  has  been  used  with  r,(490)  being  set  to  be  0.01  and  the  sun 
is  at  the  zenith. 

If  the  separation  between  the  wavelengths  selected  is  not  large  enough,  the  system  will  be  insensitive  to  the  spectral 
properties  of  the  diffuse  attenuation  coefficient,  K(X),  of  the  water  and  poor  results  will  be  obtained.  For  a  separation  of 
20nm  Fig.  5  indicates  that  480  and  500nm  are  nearly  "balanced"  wavelengths  for  K(490)  »  0.067  type  water.  The  sensitivity 
of  measurement  at  these  two  wavelengths  to  change  in  water  type,  K  (490),  is  demonstrated  in  Fig.  6.  The  value  "Q"  is  a 
"figure  of  merit”  where: 

Q  -  2.3026  •  [log  R2  -  log  R,]/[K(490)2  -  K(490),]  •  z  (10) 

and 

R  *  Ez$z)  /  Ej(Xi)  .  (11) 

The  statement  "if  AK(490)  -.001  then  AR  »  1.019"  in  Fig.  6,  indicates  that  a  change  in  K(490)  of  .001  will  result  in  a 
change  in  the  ratio  of  the  irradiances,  R  -  E,(500yE,(480),  at  depth  z  (100  meters  in  this  case)  of  1.9%.  It  can  also  be 
inferred  from  this  that  an  error  in  the  determination  of  the  ratio,  R,  of  1.9%  will  cause,  approximately,  an  error  or  .001  in  the 
computation  of  K(X).  In  a  similar  manner  but  using  wavelength  separations  of  50  nm  and  100  nm,  the  results  shown  in 
Figs.  7  &  8  are  obtained.  It  is  seen  that  increasing  the  separation  markedly  improves  the  sensitivity  to  change  in  the  water 
quality.  This  improved  performance  comes  at  the  cost  of  depth  capability.  Again  using  Fig.  5  we  can  estimate  the  depth 
limitation;  i.e.,  the  depth  at  which  the  irradiance,  E,(X),  arriving  at  one  or  both  photometers  is  0.0002p.W  ■  cm  2  •  nm*1. 
To  summarize: 

For  W«ter  Type  K(490)  -  0.067 


X](am) 

X2(nm) 

SeniiUviiy 
change  in  DR  with 

0.001  change  in  K(490) 

Approximate 
Depth  LimitaUon 
(meters) 

480 

500 

1.019 

200 

460 

510 

1.049 

180 

440 

540 

1.092 

160 

The  it  the  level  M  which  <  photometer  employing  a  high  quality  photomultiplier  tuba  having  ae  irradiance  collector  with  a  3/4  inch  diameter,  a  band  pans  of  lOnm,  and  a 
uaoanaaaioo  through  the  optical  path  of  1%,  might  be  expected  to  have  a  ii  goal  to  ooieenlioiaexcen  of  ten. 


The  above  selection  process  could  be  used  to  obtain  maximum  depth  capability  for  any  one  water  type.  Using  460  and 
SlOnm,  Fig.  9  explores  how  the  E,(460)  and  Ez(510)  values  and  their  ratio  would  change  with  depth  in  K(490)  «  0.067  type 
water.  Since  K(460)  -  0.0759  and  K(510)  -  0.0739  are  not  greatly  different,  the  ratio  does  not  change  much  with  depth. 
Figures  10  and  1 1  are  similar,  using  the  same  two  wavelengths,  but  in  water  types  K(490)  =  0.038  (Jerlov  Type  IB)  and 
K(490)  »  0.115  (Jerlov  Type  IE).  Note  that  the  arbitrary  irradiance  limit  of  0.0002pw  ■  cm2  •  nnv1  is  reached  at  a  depth  of 
about  100  meters  at  X»460nm  when  K(490)  »  0.1 15  and  has  not  been  reached  at  200  meters  when  K(490)*0.038. 

The  selection  of  wavelength  based  on  the  second  criteria,  minimum  sensitivity  to  system  error  (primarily  wavelength),  leads 
to  a  different  set  of  wavelengths  and  some  loss  in  depth  capability.  By  wavelength  error  is  meant  the  difference  between  the 
wavelength  at  which  we  think  the  radiometer  is  making  the  measurement,  (the  assigned  wavelength  used  in  the 
computations),  and  the  actual  effective  or  mean  wavelength  to  which  the  radiometer  responds.  For  instance  if  the  radiometer’s 
effective  response  was  at  492nm  and  we  used  490nm  in  the  calculations  the  wavelength  error  would  be  -2nm. 

To  acheive  minimum  sensitivity  to  wavelength  error  we  must  find  a  wavelength  where  K(X)  (and  to  a  much  lesser  extent 
E0(X),  tR(X),  t0(X)  and  ta(X))  is  not  changing  rapidly  with  wavelength  over  the  channel  band  width.  We  also  need  to  do  this  at 
two  wavelengths  which  are  far  enough  apart  and  in  regions  which  will  give  good  responsivity  to  the  quality  of  the  water. 
Observation  of  the  K(X)  vs  X  curves,  Fig.  4,  is  helpful  in  the  selection  for  this  criteria  and  also  shows  the  difficulty  in  such  a 
selection.  Certainly  the  selection  of  the  "optimum”  pair  of  wavelengths  is  highly  dependent  upon  water  type. 

To  illustrate  the  difference  in  sensitivity  to  wavelength  error  between  wavelengths  selected  with  this  in  mind  as  opposed  to 
wavelengths  selected  to  obtain  maximum  depth  capability,  the  use  of  wavelengths  420  and  490nm  will  be  compared  with  460 
and  510nm  when  the  water  type  is  again  K(490)»0.067.  In  the  figures  starting  with  Fig.  12,  the  solid  lines  are  the  "true'' 
values;  the  points  indicated  by  diamonds  (0),  are  values  which  would  be  calculated  if  the  system  were  perfect,  i.e.  had  no 
errors  including  the  values  used  for  E0(X),  tR(X)  and  x0(X);  and  crosses  (+),  are  values  which  would  be  calculated  if  the  error 
shown  in  the  heading  under  "error  factors"  is  introduced.  The  small  error  of  the  "perfect"  system  is  due  to  the  assumption, 
made  to  allow  the  inverse  problem  to  be  solved,  that  the  diffuse  transmission  properties  of  the  aerosol,  t,(X)  are  spectrally 
independent,  i.e.  x,(X)  is  the  same  at  all  wavelengths.  Two  "types"  of  aerosol  are  used:  (1)  haze,  t,(490)-0.10  and  a- 1.00, 
and  (2)  overcast,  x.(490)- 10-20  and  a»0.00.  The  type  of  aerosol  is  indicated  on  the  graphs  by  the  numbers  1,&  2  on  the  right 
hand  side.  In  all  cases  30*  is  used  for  the  solar  zenith  angle  and  lOnm  for  the  radiometric  bandpass.  Frequently  we  would  like 
to  know  the  atmospheric  and  water  attenuations  at  a  wavelength  other  than  either  of  the  two  at  which  the  measurement  is 
made.  In  these  figures  two  wavelengths  have  been  used  for  the  measurement  wavelengths,  the  problem  solved  at  one  of  these 
wavelengths,  both  with  and  without  the  applied  error,  and  then  the  results  transferred  to  a  third  "wavelength  of  interest"  X*,  via 
the  model(s).  The  figures  have  three  graphs:  (1)  transmission  through  the  total  atmosphere.  This  is  the  vertical 
transmission,  i.e.  for  unit  air  mass.  It  is  shown  for  both  atmospheric  types.  (2)  Logarithm  of  the  transmission  through  the 
total  path.  This  is  the  vertical  transmission  through  the  atmosphere,  the  air-water  interface  and  the  water  column.  It  is 
shown  for  both  atmospheric  types.  (3)  The  ratio  of  the  calculated  transmission  to  the  true  transmission  for  the  total  vertical 
path.  It  is  shown  for  one  atmospheric  type*.  All  of  these  values  are  plotted  as  a  function  of  depth,  and  are  the  results 
obtained  for  the  third  wavelength,  the  "wavelength  of  interest". 

In  Fig.  12  the  photometric  measurement  wavelengths  are  420  and  490nm  and  the  wavelength  of  interest  is  459nm.  A 
wavelength  error  of  +2nm  has  been  applied  to  the  photometric  measurement  at  Xj(490nm).  The  problem  is  solved  to  obtain 
the  results  at  490nm  and  the  values  for  X-459nm  calculated.  To  calculate  the  effect  of  the  applied  error,  values  indicated  by 
crosses  (+),  the  problem  is  solved  using  492nm,  the  erroneous  wavelength,  in  place  of  490nm,  the  true  wavelength.  Figure 
13  is  a  repeat  of  Fig.  12  except  that  the  error  applied  to  Xj  is  -2nm.  Figuresl4  and  15  show  the  effect  of  the  same  wavelength 
error  offsets  when  460nm  and  510nm,  the  wavelengths  arrived  at  previously  to  obtain  greater  depth  capability,  are  used  for  the 
measurement  wavelengths. 

Comparison  of  Fig.  12  &  13  with  Fig.  14  &  IS  demonstrates  that  the  sensitivity  to  wavelength  error  is  highly  dependent  on 
wavelength  selection.  The  penalty  for  lower  wavelength  sensitivity  is  loss  of  depth  capability.  In  this  case  the  depth  limit, 
(K  (490)«0.067  type  water  and  again  using  0.0002|iW  •  cm-2  •  nnr1  for  the  photometer's  limit  of  sensitivity),  is  about  140 
meters  at  420nm  when  420nm  and  490nm  are  used  and  about  180  meters  when  460nm  and  510nm  are  used.  This  sensitivity 
to  wavelength  error  becomes  even  more  drastic  in  more  turbid  waters.  Figures  16  through  19  is  another  case  with  the  same 
magnitude  of  wavelength  error  applied  for  water  type  K(490)«0.248  (Jerlov  Type  2).  For  maximum  depth,  527nm  and 
577nm  would  be  used  (Fig.  18  &  19)  and  the  depth  limit  would  be  about  64  meters.  For  the  other  wavelengths  used,  420  & 

•Ths  results  presented  this  wsy  an  newly  identical  for  all  atmosphere  types. 


550nm,  the  depth  limit  is  about  33  meters  (see  Figs.  16  &  17). 


4.0  OTHER  SOURCES  OF  ERROR 

4.1  RADIOMETRIC  ERROR  (other  than  wavelength) 

Two  measurements  of  the  absolute  downweliing  irradiance  at  some  depth,  z ,  are  required.  An  error  in  these  measurements 
produces  a  corresponding  error  in  the  determination  of  the  atmospheric  and  water  diffuse  attenuation  properties.  The 
radiometric  error  is  the  sum  of  the  error  in  calibration,  drift  in  sensitivity  with  time,  and  change  in  sensitivity  with 
temperature.  The  effect  of  this  error  is  a  function  of  the  two  wavelengths  used  and  in  general  is  lessened  by  increasing  the 
separation  between  the  wavelengths.  This  is  illustrated  in  Fig.  20  where  Xta420nm  and  X2*470nm,  the  water  is  pure 
seawater,  K(490)=0.022,  and  an  error  of  -5%  is  applied  to  the  irradiance  measurement  at  X^  and  in  Fig.  21  where  the  same 
error  is  applied  but  the  wavelengths  used  are  420nm  and  330nm.  Figures  22  &  23  are  similar  to  Figs.  20  &  21  but  for  water 
type  K(490)=0.125.  The  effect  of  this  type  of  error  is  independent  of  depth  and  essentially  independent  of  water  type. 

4.2  DEPTH  ERROR 

A  one  meter  offset  in  the  depth  value  used  in  the  computations  will  cause  the  errors  shown  in  Figs.  24  &  23  for  water  types 
K(490)=0.038  and  0.250  respectively.  The  errors  caused  by  a  1%  non-linearity  in  depth  are  shown  in  Fig.  26  &  27  for  the 
same  water  types.  In  all  cases  420nm  and  530nm  were  used  for  the  two  photometric  wavelengths.  An  offset  produces  errors 
in  TA,  total  atmospheric  transmisssion,  and  T,  total  transmission  atmosphere  and  water,  which  quickly  reach  an  asymptote 
and  become  essentially  independent  of  depth.  A  non-linearity,  of  course,  produces  an  error  which  is  proportional  to  depth. 

4.3  SOLAR  ZENITH  ANGLE  ERROR 

The  error  caused  by  using  an  incorrect  zenith  angle,  0„  for  the  sun  is  a  function  of  the  cosine  of  the  zenith  angle.  It  is  small 
for  high  suns,  (small  zenith  angle)  and  large  for  low  suns,  (large  zenith  angle).  Two  examples  are  given  in  Figs.  28  &  29  for 
0,a3O*  and  0,=6O*  both  with  an  applied  error  of  +5*.  This  error  is  independent  of  depth  and  water  type. 


4.4  ERROR  INTRODUCED  BY  THE  MODEL  USED 
FOR  SEA  WATER  OPTICAL  PROPERTIES 

To  use  this  system  it  is  necessary  to  use  a  relationship  which  allows  the  determination  of  the  diffuse  attenuation  coefficient 
of  the  water  at  a  second  wavelength,  KfX^),  from  knowledge  of  the  value  at  some  other  wavelength,  K(X,)  (Ref.  1  &  2).  For 
flexibility  it  is  most  desirable  that  this  can  be  done  between  any  two  wavelengths  within  a  span  of  practical  values.  The 
model  used  for  this  purpose  is  presented  in  detail  in  Ref.  2.  To  obtain  good  results  over  the  long  water  path  lengths  desired, 
it  is  necessary  to  determine  accurately  the  K(X)  values  involved.  The  prospect  of  doing  this  depends  to  a  large  extent  on  how 
well  this  model  matches  the  spectral  character  of  the  water  in  which  the  system  is  being  used.  The  data  on  which  the  model 
is  based  is  not  very  precise.  It  is  subject  to  usual  experimental  errors  and  adversely  affected  by  difficult,  at  sea,  environmental 
operating  conditions.  It  can  be  hoped  that  the  fairly  large  data  base  used  in  forming  the  model,  resulted  in  statistically 
significant  smoothing  of  the  data  base  and  that  the  model  is  satisfactory  for  this  application. 


The  correctness  of  the  model  and  the  uniformity  of  fairly  clear  deep  oceanic  water  such  that  the  model  will  match  all  such 
waters  can  only  be  determined  by  future  experimentation  and  testing.  The  model  remains  a  potential  source  of  error,  the 
magnitude  of  which,  at  this  time,  is  unknown. 


4.5  SOLAR  AND  ATMOSPHERIC  CONSTANTS 

The  sun  is  used  as  a  known  source  of  irradiance  outside  of  the  atmosphere.  The  absolute  spectral  level  of  this  irradiance, 
Eq(X),  is  a  primary  input  for  this  system.  Much  effort  has  been  put  into  determing  the  correct  values  for  this  spectral  solar 
irradiance.  It  is  believed  the  values  of  Neckel  and  Labs  (1981)  are  correct  within  approximately  one  percent  However,  this 
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must  be  considered  a  potential  source  of  error.  Figures  30  &  31  show  the  results  of  a  +5%  error  in  Eo(X)  at  420nm  using 
420nm  and  470nm  for  the  measurement  wavelengths.  Figure  30  is  for  pure  sea  water,  K(490)=0.022  and  for  Fig.  31 
K(490)=0.125.  Figures  32  &  33  are  for  the  same  water  types  and  enor  offset  but  420nm  and  530nm  are  the  wavelengths 
used.  This  type  of  error  is  independent  of  depth  and  water  type  and  has  less  effect,  at  least  on  the  atmospheric  transmission 
calculated,  when  the  wavelengths  used  are  farther  apart 


5.0  DIFFUSE  ATTENUATION  COEFFICIENT  OF  THE  WATER 

The  ability  to  determine  the  effective  diffuse  attenuation  coefficient,  K(X),  of  the  water  at  any  wavelength  of  practical  interest 
from  the  measurement  of  downwelling  irradiance  at  two  wavelengths  is  a  significant  secondary  objective  of  this  type  of 
system.  Examples  are  given,  all  using  420nm  and  530nm  for  the  measurement  wavelengths  and  the  results  shown  for  the 
computed  K(A)  at  459nm.  The  computed  K(X)  vs  depth  for  water  type  K(490)=0.067  with  an  applied  error  of  +2nm  at 
X[(420nm)  is  shown  for  two  types  of  atmosphere  in  Fig.  34.  Figure  35  has  the  same  error  applied  at  420nm  but  for  water 
types  K(490)*0.038  and  K(490)»0.152  with  atmosphere  No.  2  (overcast).  Figures  36  through  41  all  show  examples  using 
these  same  two  water  types.  In  Fig.  36,  the  error  applied  makes  the  irradiance  measurement  at  420run  low  by  5%. 
Figure  37  has  a  depth  offset  error  of +1  meter  at  all  depths  and  Fig.  38  has  a  depth  non-linearity  of  +5%.  The  effect  of  a  +5° 
error  in  the  solar  zenith  angle  is  shown  in  Fig's.  39  &  40;  with  the  "true"  zenith  angle  30°  in  Fig.  39  and  60°  in  Fig.  40. 
Figure  41  has  an  error  of  -5%  in  the  assumed  solar  irradiance  at  420nm. 

When  a  fairly  long  path  length  is  used  (2-4  attenuation  lengths;  transmission  about  2  to  10%),  the  theoretical  computations 
produce,  even  with  these  applied  errors,  values  which  are  closer  to  the  real  value  than  can  be  expected  from  the  usual  in-situ 
methods  normally  employed  to  get  this  type  of  data.  How  well  this  works  in  the  real  world  needs  to  be  determined.  The 
accuracy  of  the  K(X)  calculated  is  subject  to  degradation  by  all  the  errors  discussed  in  previous  sections.  It  is  largely  the  error 
involved  in  determining  K(X)  which  will  cause  an  error  in  the  calculated  atmosphere  and  water  path  transmissions. 


6.0  SUMMARY  AND  COMMENTS 

|  Wavelength  Selection  •  There  are  two  considerations  in  selecting  the  two  wavelengths:  (1)  maximum  depth  at  which  the 

i  photometers  will  have  adequate  signal  and  (2)  sensitivity  of  the  system  to  radiometric  and  other  error  sources.  These  two 

|  criteria  do  not  lead  to  the  same  wavelengths.  Selection  of  "optimum  wavelengths"  to  satisfy  one  has  an  adverse  effect  on  the 

i  other  and  both  are  highly  dependent  cm  the  type  of  water.  There  is  no  one  "best  choice".  The  choice  is  a  matter  of 

|  compromise,  judgement  and  pratical  considerations.  Large  separation  between  the  wavelengths  in  most  cases  will  improve 

system  performance. 

|  Radiometrv  -  If  adequate  system  performance  is  to  be  achieved,  good  radiometric  calibrations  and  stability  are  required.  The 

I  absolute  responsivity,  linearity  and  "effective  mean  wavelength"  for  both  radiometers  must  be  known.  The  radiometers  need 

to  be  stable  with  respect  to  time,  environment,  and  temperature.  The  spectral  filters  may  degrade  with  time  and  the 
sensitivity  of  the  radiometer  may  change.  This  suggests  that  periodic  "re-calibration"  should  be  done.  The  sensor, 
particularly  if  a  photomultiplier  tube,  may  suffer  a  change  in  spectral  response  with  changes  in  temperature.  Some  method 
may  be  required  to  stablize  or  correct  for  such  changes  in  sensitivity.  A  stable  internal  reference  light  source  could  be  used  to 
correct  the  readings  for  changes  in  sensitivity. 

Errors  -  The  system  can  be  very  sensitive  to  wavelength  and  the  wavelength  at  which  the  measurements  are  made  must  be 
closely  known  to  avoid  possible  large  errors.  Absolute  radiometric  error  has  less  significance  if  both  radiometers  are  in  error 
in  the  same  direction.  The  inherent  error  of  the  two  wavelength  system  is  small  relative  to  the  probable  error  caused  by  other 
|  sources  such  as  radiometric  measurement  accuracy,  the  values  used  for  the  other  inputs,  and  the  fidelity  of  the  K(X)  model 

used. 

Water  Properties  -  The  system  can  be  a  useful  tool  for  the  determination  of  the  spectral  diffuse  attenuation  properties  of 
natural  oceanic  waters. 

|  Atmosphere  •  It  appears  that  an  estimate  of  the  optical  depth  of  the  atmosphere  can  be  obtained  which  is  adequate  for 

j  assessing  the  total  atmospheric  transmittance  loss. 

i 
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Total  Transmission  •  When  the  total  spectral  diffuse  transmission,  i.e.,  the  transmission  through  the  total  atmosphere,  the  air- 
water  interface  and  the  water  path,  is  examined,  it  becomes  apparent  that  for  any  appreciable  water  path  (i.e.,  depth)  the 
attenuation  through  the  water  is  predominate  over  the  loss  in  the  atmosphere.  In  estimating  the  total  attenuation  it  appears 
that  this  type  of  system  may  produce  very  satisfactory  results. 

Final  Comment  •  The  inherent  error  in  the  computations  which  results  from  the  assumption  that  the  portion  of  the 
atmospheric  attenuation  caused  by  the  aerosol  is  constant  at  all  wavelengths  is  very  small.  For  atmospheres  which  are  not 
cloudy,  heavily  overcast,  or  foggy  (ta<l),  the  term  for  the  aerosol  in  Eq.  (1)  for  the  atmospheric  transmittance  given  in 
section  2.1,  i.e.  l/2(l-g)ra(X)  is  approximately  1 16  ■  xa(X)  and  has  a  minor  effect  on  the  results.  For  aerosols  of  greater 
densities  a  is  very  small  and  the  aerosol  attenuation  though  huge  is  essentially  independent  of  wavelength. 

All  that  has  been  presented  in  this  report  assumes  perfection,  other  than  for  the  errors  intentionally  introduced.  In  a  real 
system  the  final  results  will  be  the  summation  of  all  sources  of  error.  The  goodness  of  the  match  between  the  K(X)  water 
model  and  the  actual  waters  is  a  source  of  error  the  magnitude  of  which  is  yet  to  be  determined. 


7.0  RECOMMENDED  WAVELENGTHS 

For  a  field  test,  five  photometers  were  available  allowing  the  choice  of  five  different  wavelengths.  Two  of  those  photometers 
had  a  linear  response  with  auto-ranging.  The  remaining  three  had  a  pseudo-logarithmic  response. 

Using  the  approach  presented  in  this  report,  five  wavelengths  were  selected:  420, 470, 490, 530,  and  S50nm  with  420nm  and 
530nm  assigned  to  the  linear  photometers.  Table  A  gives  the  selected  combinations  and  lists  the  water  types  for  which  it  is 
recommended  these  wavelengths  be  used.  The  first  choice  for  wavelength  pair  is  420  and  530nm.  Table  B  gives  an  estimate 
of  certain  errors  that  might  be  anticipated  for  a  2nm  wavelength  calibration  offset,  when  using  these  two  wavelengths  in  five 
different  types  of  water.  Table  C  gives  the  second  choice  of  wavelength  pairs  which  could  be  used  to  increase  depth 
capability.  Again  the  same  estimates  of  error  are  given  for  the  five  types  of  water.  The  error  values  in  Table  C  can  be 

directly  compared  with  the  values  in  Table  B  where  Ac»459.  Tables  B  and  C  also  give  the  depth  limits  which  can  be 
anticipated  for  the  various  pairs  of  wavelengths. 

A  photometric  error  of  5%  in  the  measurement  of  the  imdiance,  EL(X),  can  cause  the  resultant  errors  shown  in  Table  D  for  the 
various  recommended  wavelength  combinations.  Similarly  the  effect  of  a  1  meter  offset  in  depth  is  shown  in  Table  E.  Table 
F  gives  the  errors  which  can  result  from  a  +5%  offset  in  solar  zenith  angle.  This  error  can  be  severe  for  large  zenith  angles 
(low  sun).  The  errors  indicated  in  Table  G  can  result  from  a  5%  error  in  the  value  used  far  E^X),  the  extra-terrestrial  solar 
irradiance  arriving  outside  the  atmosphere. 

Table  H  is  a  listing  of  the  solar  irradiance,  EofytW  •  cm*2  •  nm*>),  outside  the  atmosphere  and  the  optical  depths  Xr(X) 
(Rayleigh),  xQ(X)  (ozone)  and  Xq(X)  (aerosol);  calculated  as  discussed  in  this  report  A  lOnm  bandpass,  uniform  response 
from  X-Snm  to  X+5nm  and  a  value  of  xa(490>0.01  (clear  atmosphere)  was  used. 
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The  depth  nnge  given  is  for  the  conespooding  water  type  range  and  is  the  depth  at 
which  the  inadiance,  Ej(X),  a  either  or  Xj  Tint  teaches  0.0002|aW  car1  oar'  with 
a  dear  atmosphere  and  the  tun  a  the  zenith. 


Tnbln  B.  Wavelength  error,  Xj  -420,  X2-530nm.  Maximum  error  caused  by  a  ±  2nm  offset  in 

wavelength  calibration. 
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Notes:  (1)  Depth  a  which  Et(X)-0.0002  |iW  cm1  -  am’1  (dev  atmosphere,  sun  a  the  zenith). 

(2)  Transmittance  through  atmosphere  a  wavelength  Xe 

(3)  Transmittance  through  total  path  (atmosphere  ♦  water)  a  wavelength  \ 

(4)  Wavelength  and  direction  (±)  a  which  2nm  offset  causes  maximum  error 


Table  A.  Recommended  wavelengths. 


Xi  (am) 

Xt  (nm) 

Water  Type,  K(490),  Range  (nr1) 

Depth  Range(m) 
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260  —  33 

420 
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.022  —  .046 

620  —  229 
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197  —  97 
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550 

.125  —  250 

109  —  55 

NOTES:  The  combination  X,«420  and  Xj-530  is  the  Tint  reccommeaded  choice.  The 
other  combinations  should  be  used  only  if  necessary  to  obtain  adequate  signal. 


-II. 


V ■  n  *' 


r  v 


V  •  f* 


I, 


Table  C.  Wavelength  error  for  X,  and  Xj  in  table. 

Maximum  error  caused  by  a  2nm  offset  in  wavelength  calibration  using  X,  and  Xj  to  obtain  greater  depth  capability. 
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Nous:  (1)  Depth  at  which  E,(X)«0.0002  |iW  -cm'1  -  nm'1  (clear  atmosphere,  sun  at  the  zenith). 

(2)  Transmittance  through  atmosphere  at  X(a459nm. 

(3)  Transmittance  through  total  path  (atmosphere  ♦  water)  at  Xt»459nm. 

(4)  Wavelength  end  direction  (±)  at  which  2nm  offset  causes  maximum  error. 


Table  D.  Maximum  error  caused  by  a  5*  offset  in  the  measured  oiadiaoce.  E^X),  at  Xt  or  X^ 

Calculated  at  Xc-459nm 
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(nm)  (nm)  from  to  in  !%,(%)  in  T*.  (%) 


Table  E.  Error  caused  by  a  1  meter  offset  in  depth.  Calculated  at  \.-459nm. 


0 

(am) 

— 

**  ! 
(nm)  i 

Water  Type  K(490) 
nr1 

■ 

Enor  in  j 

(%> 

0j 

Error  in  T  j 
(*>  | 

420 

— 

530 

.022 

8  i 

3  i 

420 

530 

.046 

8 

7 

I  420 

530 

.077 

8 

11 

1  420 

530 

.125 

8 

18 

Ea 

250 

8 

39 

420 

490 

.022 

3 

2 

420 

490 

.046 

3 

5 

470 

530 

.046 

11 

470 

.125 

11 

16 

490 

550 

.125 

13 

13 

490 

550 

250 

13 

33 

Notes:  (1)  X,,  X,  are  the  two  wavelengths  at  which  the  inadiance,  E,  (X),  are  measured 
(The  man  are  calculated  for  results  transferred  to  Xc-4S9nm). 

(2)  Transmittance  through  atmosphere  at  Xt^4S9nm. 

(3)  Transmittance  through  total  path  (atmosphere  *  water)  at  Xc-459nm. 

(2),  (3)  Error  is  independent  of  depth. 


Table  F.  Error  caused  by  a  «-5%  offset  in  solar  aenith  angle,  8, 


8,  (deg.) 

ERROR  IN  TA(%) 

Ham  Overcast 

ERROR  IN  T(%) 

Hare  Overcast 

25 

3 

10 

3 

10 

30 

5 

19 

5 

19 

45 

7 

32 

8 

33 

60 

9 

55 

11 

59 

75 

10 

133 

25 

170 

Notes:  (1)  Error  is  independent  of  depth. 

(2)  The  wavelength  at  which  the  computations  are  made,  X* 
has  negligible  effect  on  the  error  (*±0.2%). 

(3)  The  measurement  wavelengths,  \\  X*  used  have  a  small 
effect  on  the  resultant  error  (-  ±0.4%). 

(4)  The  error  is  independent  of  type  of  water. 

(5)  Ta  m  transmittance  through  the  atmosphere. 

(6)  T  ■  transmittance  through  the  total  path  (atmosphere  +  water). 


Table  G.  Error  caused  by  a  ±5%  offset  in  the  solar  inadiance,  EJX). 
Calculated  at  X,«459nm. 


Xj  (nm) 

mm 

®1 

Error  in  T'(%) 

420 

530 

9 

3 

420 

490 

14 

3 

470 

530 

13 

7 

490 

550 

14 

10 

Notes:  (1) 

(2) 

(3) 

(2).  (3) 
(3) 


Xf,  A*  are  the  two  wavelengths  at  which  the  imdiance,  E,  (X),  are  measured. 
(The  errors  are  calculated  for  results  transferred  to  X^SSnm). 

Transmittance  through  atmosphere  at  Xc«459nm. 

Transmittance  through  total  path  (atmosphere  +  water)  at  Xe<-459nm. 

Errors  in  ”FA  and  T'  are  independent  of  depth. 

Errors  in  T*  can  be  larger  or  smaller  if  another  "wavelength  of  interest",  X,., 
is  used. 


c* 


MODEL  MLC761625 
BMNOPHii  10  n» 


I 

cs 


|R 

iff 


I 

g 

g 


i 

& 


1 1  i>e  ■  «  0  01®  •*  »  i  298 


.<  r.M 

to 

fOr 

T  4 

aie 

17® 

99 

® 

1334 

0 

0022 

0 

0 1 26 

■»;3 

173 

29 

0 

1389 

0 

0826 

0 

0124 

420 

172 

62 

0 

1433 

0 

0031 

0 

0122 

423 

163 

80 

0 

1372 

0 

0O37 

8 

0120 

430 

162 

49 

0 

1309 

0 

9044 

0 

01  13 

433 

160 

47 

0 

1230 

0 

0032 

0 

0117 

440 

133 

33 

0 

!  194 

0 

0061 

0 

01  If 

443 

192 

73 

0 

1  141 

0 

0071 

0 

0113 

450 

200 

19 

0 

1091 

0 

0063 

0 

01  12 

433 

203 

39 

0 

1044 

0 

0097 

0 

0110 

460 

203 

37 

0 

0999 

0 

0112 

8 

0109 

463 

200 

82 

0 

0937 

8 

0130 

0 

0107 

470 

199 

14 

0 

0917 

0 

0130 

0 

0186 

473 

201 

19 

0 

0879 

0 

0173 

0 

0104 

400 

201 

86 

0 

8343 

0 

0198 

0 

0l®3 

405 

194 

04 

0 

0809 

0 

0227 

0 

0101 

490 

139 

03 

0 

0776 

0 

0258 

0 

0100 

495 

193 

63 

0 

0743 

0 

0293 

0 

0099 

300 

192 

48 

0 

0716 

0 

0332 

0 

0097 

303 

191 

29 

0 

0633 

0 

0374 

0 

0O96 

510 

191 

43 

9 

0661 

0 

0421 

0 

0893 

313 

132 

72 

0 

0636 

0 

0471 

8 

0094 

320 

131 

62 

0 

0612 

0 

8325 

0 

0093 

323 

136 

36 

0 

0589 

0 

8383 

0 

0O91 

33® 

1 S7 

38 

8 

0367 

0 

0644 

0 

0096 

335 

193 

43 

0 

0346 

0 

8709 

0 

0089 

540 

183 

90 

0 

0326 

0 

0777 

0 

0088 

543 

183 

16 

0 

0307 

0 

0847 

0 

0087 

33« 

136 

29 

0 

0489 

0 

0919 

0 

0086 

533 

133 

63 

0 

0472 

0 

0992 

0 

0685 

340 

182 

33 

0 

0453 

0 

1065 

0 

0084 

363 

173 

42 

0 

0439 

0 

1136 

0 

0083 

370 

174 

12 

0 

0424 

0 

1204 

0 

0082 

573 

191 

7* 

0 

0409 

0 

1266 

0 

0881 

300 

133 

24 

0 

0395 

0 

1321 

0 

0080 

Fig.  lb.  Table  showing  the  effective  extra-terrestrial  solar  irradiance  is  a  lOnm  band  centered  at  the  wavelength  in  column  1, 
and  the  optical  depths  tR(X)  (Rayleigh  or  molecular  component),  t0(X)  (oione),  and  t,(X)  (aerosol).  Uoiform  sensor  response 
over  the  pass  band  and  ®  t,(490)  -  0.01  were  assumed. 
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SDEV 

.8026 

Fig.  2.  Ozone  spectral  absorption.  Data  and  model. 
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Fig.  4.  Optimum  wavelengths  for  maximum  depth  vs  K  (490). 
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Fig.  S.  Depth  limit  (for  0.2}iw  •  cm'2  *  tun-1)  vs  wavelength  when  K(490)-0.067m't 
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Fig.  6.  Sensitivity  of  proposed  method  to  changes  in  water  type  (K(490))  when  wavelengths  are  Xj«480  and  X2-S00nm. 
Depth  in  example  is  100m.  For  definition  of  Q  see  text. 
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Fig.  7.  Sac bb  aa  Fig.  6  except  wavelength*  are  X,>460  and  Xj-SlOnm  (AX«50nm). 
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Flf.  8.  Same  u  Fig's.  6  &  7  except  wavelengths  are  X,«440  and  X2*540ntn.  (AX  »  lOOnm). 
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Fig.  9.  Sensitivity  to  changes  in  depth  with  wavelengths  Xj-460  and  X}«510na.  Water  type  K(490)>0.067. 
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Flf.  10.  Sum  u  Fig.  9  except  K(490)-0.038. 
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Flf.  11.  Sidm  u  Fig's.  9  ft  10  except  K(490)-0.115. 
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Fig.  12.  Sensitivity  to  wavelength  error.  K(490)«0.067,  f,  -  30*.  Curve  tabled  "1“  for  haze,  t^O.IO,  X-1.00.  Curve 
labeled  "2"  for  ovcrcait  t, -10-20,  X=0.0.  Results  shown  for  X=459nm  and  unil  air  mass.  Measurement  wavelengths  X,«420 
Xj-490nm.  Error  >2nm  in  X,- 
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Fig.  IF.  Same 


Fig.  17.  Suae  u  Fig.  16  except  eiror  it  -2am  in  X 


Fig.  16  except  X,«527nm,  X2-577nm.  Error  is  +2nm  in  X 


Fig.  19.  Same  as  Fig.  18  except  error  is  -2am  ia  V 


Fig.  20.  Sensitivity  to  Radiometric  error,  K(490)>0.022m  ,l  A.,.420,  X2-470nm.  Error  is  -5%  in  E^X,).  Other  conditions 

as  in  Fig.  12. 
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24.  Seatitivity  to  errors  in  depth  of  u/w  sensor.  K(490)»0.038m'',  X|»420,  Xj*530.  Error  is  1  meter  offset  in 
ir  conditions  is  in  Fig.  12. 
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fig.  34.  Sensitivity  of  K  determination  to  errors  in  wavelength,  K(490)«0.067m',l  *.,-420,  X2-530nm.  Error:  +2nm  in  X, 
(a)  Haze  (t,-0.1,  a- 1.0),  (b)  Overcut  (t,- 10-20,  a-0).  Other  conditions  u  in  Fig.  12. 
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Fig.  35.  Same  as  Fig.  34  except;  (a)  K(490)-0.038m‘.  (b)  K(490)-0.152m  ‘. 
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Fig.  36.  Sensitivity  of  K  determination  to  radiometric  error.  X,«.420,  X2-530nm.  Error:  -3%  in  Et(X,).  (a)  K(490)=0.038m 
(b)  K(490)«0.152nr1.  Other  conditions  at  in  Fig.  12. 
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Fig.  38.  Same  as  Fig.  37  except  +5%  error  io  Z. 
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Fig,  41.  Sensitivity  of  K  determinations  to  errors  in  solar  spectral  inadiance.  X,»420,  X2-530.  Error:  -5%  in  E0(X,). 
(a)  K(490)-0.038m*i,  (b)  K(490)-0.1S2m'1.  Other  conditions  as  in  Fig.  12. 
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14  February  LABI 


EN-001-3  2t 


TECHN I  CAL  MEMO  RANDU.M 


SUBJECT:  Determination  of  the  Optical  Properties  of  the  Atmosphere  above 

the  Ocean's  Surface  from  Spectral  Measurements  cf  the  Natural 


Light  made  below  the  Surface. 


1 . 0  INTRODUCTION 


This  brief  study  was  made  in  response  to  the  question:  From  optical 
measurements  made  at  an  arbritary  depth  (z)  in  the  open  oceans,  can  the 
optical  character  of  the  atmosphere  existing  at  the  time  of  the  measurement (s) 


be  determined?  That  is,  can  we  conceive  of  a  technique  which,  using  the 


ambient  light  at  some  depth  in  the  ocean,  will  provide  information  from  which 
we  can  determine  the  type  of  atmosphere  (clear,  hazy,  overcast,  etc.)  existing 
overhead?  It  is  assumed  that  the  depth  (z) ,  at  which  the  measurement  is  made, 
and  the  zenith  angle  of  the  sun  (9s)  (i.e.  time  of  day  and  latitude  etc.)  are 


known. 


2 . 0  CONCEPT 


Prior  work  indicates  that  there  exists  an  apparent  spectral  relation¬ 


ship  between  the  diffuse  attenuation  coefficient,  K(\),and  the  various  "types" 


of  natural  waters  found  in  the  open  oceans.  The  attenuation,  K(\),  varies  with 
wavelength  and  its  spectral  form  is  quite  different  for  clear  water  than  for 
more  turbid  water.  If  this  relationship  is  consistent  and  known  for  the 
waters  of  interest,  then  by  making  measurements  of  the  ambient  spectral 


irradiance  at  deDth  we  should  be  able  to  determine  the  attenuation  characteris¬ 


tics  of  the  water  column.  Assuming  we  can  determine  the  diffuse  attenuation 


of  the  water  above  the  instrument  and  that  the  instrument  measures  the  absolute 


downwelling  irradiance,  the  downwelling  irradiance  at  the  surface  can  be  com¬ 


puted.  From  the  irradiance  at  the  surface  the  attenuation  through  the  atmos¬ 


phere  can  be  found  using  established  values  of  the  solar  irradiance  outside 


the  atmosohere. 


J.D  TWO  WAVELENGTH  MODEL 


A  two  wavelength  model  Is  used  to  illustrate  the  procedure.  The 
wavelengths  (')  used,  \  =  490  nm  and  '  =  520  run,  were  chosen  to  take  advan¬ 
tage  or  a  fairly  large  body  of  data  which  has  been  processed  by  the  Visibility 
Laboratory  to  yield  the  spectral  relationship  of  the  diffuse  attenuation  co¬ 
efficient  *<;V  at  these  specific  wavelengths.  These  may  not  be  the  optimum 
wavelengths,  but  they  serve  well  to  illustrate  the  mechanism  of  the  concept. 


3.1  K(M  Relationships 

A  study  at  the  Visibility  Laboratory  ^  has  developed  the 
following  relationships: 

-1.491 

+  0.022  (1) 

-1.398 

+  0.044  (2) 

Where  Lss ( ' )  is  the  subsurface  upwelling  radiance  at 
wavelength  X. 


K (490)  =  0.0883 


K  ( 520 )  =  0.0663 


Lss (443) 


Lss (550) 


Lss (443) 


Lss (550) 


L 


From  this  it  follows  that 

f 


1  0.9376 


K(520)  =  0.0663  |  11.325  x  K(490)  -  0.2492 J  + 

3.2  Propagation  of  the  Light  Downward 
Let: 

=  the  transmittance  through  the  atmosphere 
Tc  =  the  transmittance  through  the  water  surface 


0.044 


(3) 


(1)  Ref.  R.W.  Austin;  Remote  Sensing  of  the  Diffuse  Attenuation  Coefficient  of 
Ocean  Water,  Visibility  Laboratory,  29th  Symposium  of  the  AGARO  Electromagnetic 
'Wave  Propagation  Panel;  Monterey,  Calif.,  Aoril  1981. 


T  =  the  transmittance  through  the  water  to  deDth  z 

w 


jo  =  cosine  of  solar  zenith  angle 


All  trar.snittances  are  for  diffuse  light. 


Then  the  irradiance  at  depth  z  is 


EZ(.\)  =  40  *  Ta  !  !  *  TS  (  ’  TW  ( *  Eo  ' '  )  where 


the  solar  irradiance  outside  the  atmosohere. 

T  ,  .  .*£  jjT»«*  ♦  To(X,  ♦  V.\), 


Where  T  ,  T  and  T  __  are  the  optical  depths  of  the  atmosphere 

R  O  cl 

for  unit  air  mass  for  Rayleigh  (air  molecules \  for  ozone,  and  for 
the  aerosol,  respectively.  The  values  for  Tr(a)  and  T  (\)  are  fairly 
well  known  and  are  essentially  fixed,  appearing  to  vary  slightly  with 
season  and  latitude.  The  aerosol  component  has  a  wavelength  dependency 
which  appears  to  follow  a  power  law  so  that 


hr\  - 


T ( X 1 )  =  •  t(X2)  . 


The  exponent,  <*,  is  known  as  the  Angstrom  exponent  and  is  commonly 
taken  to  be  in  the  range  of  1.1  -1.4  for  fairly  clear  atmospheres  and 
it  will,  we  believe,  be  around  zero  for  very  heavy  (fog,  overcast) 
aerosols. 

Ts  =  0.98  (7) 


-K(\)  •  z 


Tw  =  e 


And  we  get 


!  (,\)  =  uo  •  (0.98) 


.  (V'>  *  V">  *  V1)  . 


E  (9) 


3.3  Solving  the  problem 


The  above  (section  3.2)  is  used  only  to  calculate  the  irradiance 
E;  which  would  arrive  at  depth  z  for  various  types  of  atmospheric 

9 

conditions  and  sea  water, to  use  for  inputs  in  solving  the  problem. 

At  this  juncture  we  are  at  depth  z  with  an  instrument  which  measures 
E^ ( ■ )  at  two  wavelengths  490  and  520  ma.  From  these  two  measurements 
we  wish  to  determine  the  type  of  atmosphere  overhead. 


Let  be  the  first  wave length, j\^  =  490  nmj  and  \0  be  the  second 
520  nmj.  From  ( 9)  we  can  get  the  ratio 


Ed)  -  —  't  (1)  -  T  (2)  +  T  (1) 
Z  UO  \  R  R  0 

-  =  e  ' 

E  (2) 
z 


-  —  T,  ( 1 )  — T  (2) 
uo  a  a 


e-  K  (1)-K  (2)|  •  Z. 


Tr  (A)  and  Tq  (A)  are  very  nearly  constant  so  let 
.  To"  (tr  (1)  -  V2'  +  Toll>  -  To<-’>) 


C1  • « 


T  d)  and  T  (2),  the  optical  thickness  due  to  aerosols  and  clouds 
are  unknowns  which  we  are  seeking  to  establish.  However,  since  A^  and 
A ^  are  close  together  and  the  Angstrom  exponent  in  Eq.  (6)  is  small 
(< 1 . 4 )  we  assume  that  they  are  equal  and  then 


( 1 )  — *  T 

a 


* 

1 . 0  . 


thus 


E  (1)  -  (k  ( 1 )  -K  ( 2 )}  •  z 

-  ~  C  •  e 

E  (2)  1 

z 

*  See  section  4  for  a  discussion  of  the  manr.itude  of  error  caused  by 
this  assumption. 
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k  = 


.  ,  e  (i; 

1.1  2 

~  z  C  *  "  < ~ ) 

1  z  ■" 


Kd)-K>2;  = 


-  In  C, 
z  1 


E  '.2) 

z 


Substituting  Eq. (3)  into  (12)we  get 


K  ;  2) “0 . 0663  (11.325) *K  (1J-0.2492 


E  (2) 
z 


0.3376 


-0.044 


=  -  lr.  C 


1  E(l) 

z 


Where  K  ( 1 ) is  a  calculated  value  as  opposed  to  X(l)  the  "true” 
value. 

From  Eq.  (13)  by  successive  approximation  K(i)  can  be  calcu- 

I 

lated  and  if  desired  K(l)  can  be  calculated  using  (3). 

» 

The  irradiance  at  the  surface  can  now  be  estimated  using  K ( A ) , 


,  1  K'U)  •  z  „  ,,, 

E’s(X)  '  05  ’  e  • 


The  calculated  atmosoheric  transmittance  is  then 


.  VA) 

TA(X)  uo  ’  E  (A) 
o 


and  since 


2  (■'■)  =  e 

A 


(A)  +  -  (A)  + 


T  (.'.) 

a 


-to  •  In  T  (a)-  t  (  ■.)  -  T  ('' ) 


=  -uo  •  In  t:  Esl  '  -  t.C)  T  (  • ) 

uo  FTT  R  -  ° 

o 
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)  =  -  UC 


In 


9 . 98 • u : 


•  e 


K  ( '  )  *  z 


(  '  ) 


Eo('> 


,(  ' 


, !  '  ) 


r 

t  * 

■  (X)  =  —  uo  •  ‘K  (A)  •  z  +  In 

CL. 

l 

Ez(X)n 

- 

ft 

0.98*  to 

E0(X) 

s 

L 

- 

The  magnitude  of  t  (X)  is  an  indicator  of  the  aerosol  density. 

cl 


For  examde: 


T  (X)  >  .1 
a 

~  .5 
>  1 


Clear  -*•  fairly  clear 
Haze  -  light  overcast 
Overcast  -*■  heavy  overcast 


4.0  COMMENTS 


4.1  The  duffuse  attenuation  coefficents,  K(X),  are  the  mean  or 
"effective"  values  over  the  total  water  column  above  depth  z.  The 
models  used  for  K(490)  and  K(520)  were  developed  using  data  roughly 
over  one  attenuation  length  (z  =r/K(X))from  the  surface  downward. 

It  has  been  assumed  that  these  relationships  hold  for  the  mean  K(\) 
over  greater  depths  (  1  -►  200  meters) .  This  needs  to  be  verified. 

4.2  The  solution  determines  the  optical  thickness  normalized  to 
one  air  mass  but  from  conditions  mainly  in  the  direction  of  the  sun 
not  directly  overhead. 

4.3  Assuming  a  non-spectrally  dependent  aerosol  in  going  from  (10)  to 
(12) (Section  3. 3) results  in  an  error,  at  this  point,  of  approximately 
1%  for  Hi  sun  and  clear  air;  2%  for  Hi  sun  and  heavy  overcast.  The 
worst  case  would  be  with  a  low  sun  and  hazy  sky  where  the  error  (at 
this  point)  for  a  sun  zenith  angle  of  80°  would  be  approximately 
10%.  A  choice  other  than  zero  might  be  slightly  better  here.  Also, 

an  itterative  process  could  be  used  where  the  "  ■  "  exponent  is  estimated 
from  the  T  found  and  the  computations  repeated, etc. 

cl 


4.4  Errors  of  ±  in  the  measurement  of  the  irradiances  E 

2(1) 

and  E  do  not  seriously  degrade  the  ability  to  estimate  the  tyw* 
z  (2) 

of  atmosphere  present. 

4.5  The  precision  with  which  the  attenuation  coefficients  K(\) 
are  determined  is  increased  with  increased  path  length  (z) .  On 
the  other  hand  over  long  path  lengths  the  effect  of  error  in  K(X) 
becomes  more  detrimental  to  the  final  results.  These  two  things 
are  somewhat  compensating. 

4.6  The  place  where  this  model  is  weakest  is  probably  in  the 

K ( A )  relationships  used.  Although  it  is  based  on  a  fairly  large 
body  of  data  (88  points  for  X  =  490  and  X  *  520)  the  precision 
and  variability  of  the  data  base,  as  is  always  the  case  with  this 
type  of  in-situ  measurement,  is  not  as  good  as  we  would  like.  We 
are  left  with  questions.  Will  this  K(X)  model  or  some  variation  of 
it  work  adequately  over  the  distances  contemplated  and  how  well 
will  it  work  when  used  under  actual  in-situ  conditions.  These 
and  other  questions  are  best,  and  perhaps  can  only,  be  answered 
by  field  work  in  the  ocean. 

4.7  The  wavelengths  used  here  are  not  necessarily  the  optimum 
wavelengths.  In  choosing  wavelengths  we  wish  to  work  with  those 
which  will  provide  the  best  determination  of  the  attenuation  of 
the  water  so  that  our  extrapolation  back  to  the  surface  is  most 
accurate.  We  must  also  not  forget  the  approximation  made  in  Section 
3.3  which  becomes  less  significant  with  less  separation  between  the 
wavelengths  used. 

4.8  Examples  of  the  two  wavelengths  scheme  have  been  worked  ud  for 
a  variety  of  input  conditions.  The  spectral  irradiances  at  depths 
z,  i.e.  E  (A  )  and  E  (X?)  ,  were  computed  using  Eq. (10)  for  the  set  of 
assumed  input  parameters  (i.e.  all  water,  atmospheric  and  geometric 
variables  are  known  and  the  underwater  irradiances  computed.  Using 

the  values  of  E  (X  )  and  E  (X  )  that  were  computed  above  and  a  knowledge 


of  depth  beneath  the  surface,  time,  date,  latitude,  and  longitude, 

aDproDriate  values  of  UC ,  T  f  A  )  ,  T  (A.),  T  (A,),  and  T  (A_)  may  be 

HI  K  2.  Ol  o  2. 

determined, and  from  Eq. (11) ,  K '  ( A  )  ,  the  estimate  of  the  diffuse 

^  I 

attenuation  coefficient  at  A^  may  be  determined.  K  (l^)  may  then 
be  inserted  in  Eq. (14)  to  obtain  an  estimate  of  the  aerosol  optical 


thickness  T  (A).  With  only  two  wavelengths,  hence  two  equations,  it 
a 


is  not  possible  to  solve  completely  for  K,  T  ,  and  *,  the  three 

a 


unknown  variables  in  the  comolete  radiative  transfer  eauation.  Hence 


the  necessitv  to  assume  T  (A  )  ~  T  (A.),  i.e.  that  the  Angstrom 

a  2  a  1  . 


exponent,  <*.  is  sufficiently  small  and  yi-  is  sufficiently  close  to 


1  that  1  -r^ 


•*  1.  This  assumption  is  not  correct  in  the  general 

case  and  in' order  to  obtain  a  solution  when  x  (A  )  x  (A  )  we  must 

a  l  v  a  2. 

resort  to  measurements  at  a  third  wavelength  for  which  we  must  obtain 
a  relationship  similar  to  that  of  Eos. (1)  and  (2).  We  would  then  have 
three  equations  of  the  form  of  Eq.  (9)  and  three  unknowns  and  could 


solve  for  X  (A  )  the  optical  thickness  of  the  aerosol  component  of 
a  l 


the  atmosphere  at  the  wavelength  of  interest  with  K(A)  and  ®  as 
byproducts . 
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Abstract 

A  study  of  the  spectral  nature  of  the  diffuse  attenuation  coefficient  of  light,  A  (A),  for  various  types  of  oceanic  waters 
has  been  performed.  These  attenuation  spectra  were  computed  from  downwetling  spectral  irradiance  data.  Ed(k).  obtained 
by  U  S.,  French  and  Japanese  investigators,  working  in  widely  separated  oceanic  regions  and  using  different  measuring 
techniques  and  equipment.  Attenuation  properties  were  calculated  over  the  spectral  region  from  365  to  700 nm  and  for 
depths  from  near-surface  to  in  excess  of  100  meters. 

Examining  the  A  (A)  data,  we  find  strong,  simple,  and  useful  relationships  exist  between  the  value  of  A  at  some 
selected  reference  wavelength.  k0 ,  and  the  value  of  A  at  some  other  wavelength  such  that 
A<X)-  .V/ (\ )  [A  (K0 )  -  A,  (A<,))  +  A.  (X)  ,  where  A.  is  the  attenuation  coefficient  for  pure  sea  water  For  oceanic 
waters  (for  example.  Jerlov  types  I  through  III)  the  relationships  are  linear.  These  relationships  appear  to  be  useful 
throughout  the  entire  spectral  range  examined  and  are  particularly  good  between  420  and  say  580  nm. 

The  significance  of  the  existence  of  such  relationships  is  that  they  allow  the  inference  of  the  spectral  attenuation 
coefficient  at  all  wavelengths  from  the  attenuation  value  at  a  single  wavelength,  and  provide  analytical  expressions  for 
modeling  the  spectral  nature  of  the  attenuation  in  ocean  and  clear  coastal  water. 

Introduction 

The  need  for  information  with  respect  to  the  optical  properties  of  ocean  water  greatly  exceeds  its  availability.  For 
example,  little  data  exists  on  the  spectral  diffuse  attenuation  coefficient  of  sea  water,  if  one  insists  on  data  based  on  in  situ 
measurements.  If  we  examine  such  attenuation  spectra,  however,  we  find  that  they  behave  in  a  well  ordered  and  regular 
fashion  over  most  of  the  spectrum,  particularly  if  we  restrict  our  attention  to  the  open  ocean  where  the  water  may  be 
classified  as  Jerlov  Types  I  through  III  (Jerlov  1976).  The  concept  of  the  Jerlov  water  types  itself  is  based  on  a 
specification  of  attenuation  spectra.  For  many  applications  the  Jerlov  classification  may.  in  fact,  provide  an  adequate 
description  of  the  spectral  nature  of  the  water  attenuation. 

We  have  reexamined  the  problem  using  new  spectral  data  as  obtained  by  a  variety  of  investigators,  working  in  vari¬ 
ous  geographical  regions,  using  different  types  of  equipment  and  varied  methods.  The  results  which  we  have  obtained  are 
similar  in  many  respects  to  those  presented  by  Jerlov  but  show  some  important  differences  in  detail.  The  method  of 
presentation  of  the  results  allows  the  user  to  estimate  attenuation  spectra  for  ocean  water  knowing  only  the  value  of  the 
diffuse  attenuation  coefficient.  A.  at  a  single  wavelength.  Thus,  the  value  of  A  at  any  wavelength  that  may  be  available 
for  a  particular  water  may  be  used  to: 

(a)  estimate  the  value  of  A  at  any  other  wavelength  of  interest. 

(b)  estimate  the  complete  spectral  characteristics  of  the  diffuse  attenuation  coefficient  for  the  water,  or 

(c)  determine  the  Jerlov  water  type. 


Discussion 

It  is  helpful  conceptually  to  partition  the  diffuse  attenuation  coefficient  into  components  associated  with  known  con¬ 
stituent  components  of  the  water.  For  example,  subtracting  A. .  the  A  associated  with  pure  sea  water,  leaves  the  contri¬ 
bution  to  A  resulting  from  all  material  suspended  and  dissolved  m  the  water,  denoted  here  as  A’  Thus  if  A  (A)  is  the 
total  diffuse  attenuation  coefficient  at  wavelength  A.  then 

A'(a)-  A<A)—  A.  (A)  (1) 


A  (A)  may  be  further  subdivided  into  contributions  from  components  due  to  say.  (a)  the  presence  of  phytoplankton 
and  the  detrital  material  that  is  associated  with  it;  (b)  to  the  absorption  of  dissolved  yellow  substance  which  is  found  in 
varying  concentrations  in  ocean  water;  and  (c)  to  the  suspended  particulates  which  do  not  co-vary  with  phytoplankton  or  its 
related  pigment  concentrations. 

Smith  and  Baker  (1978),  for  example,  used  a  model  wherein  A'  was  composed  of  two  components,  i.e.  one  propor¬ 
tional  to  concentration  of  chlorophyll-like  pigments  and  the  other.  A,  (A). "...  a  variable  representing  the  average  contribu¬ 
tion  to  spectra)  attenuation  not  directly  attributable  to  chlorophyll-like  pigment?  They  successfully  reproduced  the  general 
shape  of  measured  spectral  attenuation  spectra  using  empirically  derived  values  of  A. (A).  A, (A)  and  the  specific  spectral 
attenuation  coefficient  for  chlorophyll,  M A),  together  with  a  knowledge  of  the  pigment  concentrations. 
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If  one  attempts  to  model  the  attenuation  K'(k)  using  varying  amounts  of  two  or  three  component  spectra,  each 
dependent  on  a  concentration,  then  one  must  know  the  value  of  A'<X)  at  a  minimum  of  two  or  three  wavelengths  in  order 
to  solve  for  the  various  concentrations  and  specify  the  complete  K  (X)  spectrum.  Such  a  technique  has  the  potential  advan¬ 
tage  that  the  concentrations  of  these  constituents  may  be  determined  provided  that  a  unique  set  of  specific  attenuation 
spectra  has  been  determined.  We  have  elected  in  this  study,  however,  to  attempt  to  determine  a  single  spectral  component 
A.  "lx)  which  when  added  to  Kw  (X).  the  spectral  attenuation  of  pure  sea  water,  will  provide  the  complete  diffuse  attenuation 
coefficient  specification  of  the  water,  K  (X). 


The  Data  Base 

Downwelling  spectral  irradiance  data  from  various  sources  were  first  examined  for  general  suitability,  spectral  range, 
and  availability  of  contemporaneous  data  at  other  depths.  Irradiance  data  from  148  stations  each  having  measurements  at 
from  2  to  12  depths  were  interpolated  to  provide  values  of  irradiance.  £,  (X ) ,  at  uniform  Snm  steps.  Various  pairs  of  these 
spectra  at  each  station  were  then  used  to  compute  the  K  (X)  for  the  selected  depth  interval  using  the  relationship. 


K  (X ) 


1  .  EAk.Zi) 

-'2- -I  &(X.--|) 


(2) 


These  K  <x )  spectra  were  then  examined  for  artifacts  such  as  might  be  attributable  to  changes  in  environmental  con¬ 
ditions  during  the  acquisition  of  (he  irradiance  data,  ship  shadows,  etc.,  and  a  selection  made  of  the  spectra  to  be  used  as  a 
data  base.  No  more  than  one  spectra  was  selected  from  a  station*.  The  value  of  K  at  490  nm  was  used  as  a  classification 
index  for  the  spectra.  No  data  were  used  if  the  K (490)  value  exceeded  0.25(m-1).  This  limit  excluded  relatively  few  sta¬ 
tions  in  areas  such  as  enclosed  bays  and  in  coastal  areas  of  very  high  productivity  or  regions  of  high  terrigenous  input. 
Those  excluded  data  do  not  affect  the  generality  of  the  results  for  applications  to  the  open  ocean  and  many  coastal  regions. 
For  the  special  cases  where  K  (490)  exceeds  0.25(w~ 1 )  the  shapes  of  the  spectra  are  more  disparate  and  there  is  real  merit 
to  using  a  multicomponent  model  as  suggested  above  to  obtain  a  better  fit  to  the  data  and  to  obtain  a  measure  of  the  con¬ 
centration  of  the  several  constituents  of  the  water**. 

Table  l  lists  the  sources  of  the  data  used  in  the  present  study.  The  three  groups  of  investigators  that  acquired  the 
irradiance  spectra  were  Smith  and  his  colleagues  at  the  Visibility  Laboratory  of  the  Scripps  Institution  of  Oceanography; 
Morel  and  his  colleagues  at  the  Laboratoire  de  Physique  et  Chemie,  Villefranche-sur-Mer:  and  Okami.  Kishino  and  Sug- 
ihara  at  the  Physical  Oceanography  Laboratory  of  the  Institute  of  Physical  and  Chemical  Research  in  Japan.  The  geograph¬ 
ical  regions  include  open  ocean,  upwelling  areas,  coastal  waters  and  marginal  seas.  The  irradiance  spectra  were  obtained  at 
depths  from  I  to  118  meters.  The  K's  were  calculated  using  irradiance  data  separated  in  depth  by  from  5  to  78  meters. 
Because  the  absorption  of  water  becomes  particularly  significant  at  the  longer  wavelengths,  the  spectral  range  varied  with 
the  depth  of  the  deeper  of  the  two  irradiance  spectra  used  to  calculate  AT  (X ).  The  values  of  the  classification  index, 
K  (490),  varied  from  a  very  clear  0.025m-1  in  the  Sargasso  Sea  to  a  turbid  0.245  near  the  mouth  of  the  Mississippi  River. 
A  total  of  76  K  (x)  spectra  were  selected  for  the  data  base  from  the  148  stations  examined. 

Diffuse  attenuation  coefficient  of  pure  sea  water 

The  process  of  isolating  the  attenuation  effects  of  the  material  suspended  or  dissolved  in  the  sea  water  required  a 
knowledge  of  the  attenuation  due  to  the  base  sea  water  to  which  these  materials  were  added.  Excellent  reviews  of  the 
knowledge  of  the  attenuation  properties  of  clear  natural  waters  have  been  provided  by  Morel  (1974)  and  by  Smith  and 
Baker  (1981)  Briefly  K.  for  pure  water  may  be  approximated  by  the  sum  of  the  absorption  coefficient,  a.,  and  the  back- 
scattering  coefficient  b'n.  In  "pure"  sea  water  the  scattering  would  be  molecular  or  “Rayleigh”  with  a  phase  function  having 
equal  forward  and  backward  lobes.  The  backscattering  coefficient  was  thus  assumed  to  be  one  half  the  total  scattering 
coefficient  and  K,(k).  therefor,  determined  as 


Af.  (X )  -  a.  +  '/:  b,w  . 


(3) 


Using  a.  (x)  as  given  by  Morel  and  Prieur  (1977)  and  values  for  b,w  (X),  the  sea  water  scattering  coefficient,  presented  by 
Morel  (1974),  we  have  arrived  at  the  ff.(X)  values  listed  in  Table  2  in  the  column  labeled  “Morer.  These  values  have 
been  smoothed  and  interpolated  to  orovide  estimates  of  K,  (X)  every  Snm  (only  lOnm  intervals  are  presented  in  the 
table)  The  entries  differ  by  a  small  amount  from  those  given  by  Smith  and  Baker  in  their  Table  I  (usually  less  than  one 
unit  m  the  third  decimal  place)  due  to  procedural  differences. 

Because  subsequent  processing  of  the  data  was  to  be  performed  by  a  computer,  it  was  deemed  desirable  to  find  an 
analytic  fit  to  the  (x)  curve.  As  can  be  seen  in  Fig.  1(a),  the  curve  does  not  lend  itself  to  fitting  with  a  single  function. 


•  On  the  S C O  R  DISCOVERER  npodition  Mem  an 4  Smith  both  meaeurod  leectnl  irradiance  unfit  dilTertm  equipment  and  method!  Of  the  IT  nation  (electee  from  Moral'!  *ore 
and  the  T  iianom  Mtecied  from  Smith'!.  3  een  the  tame  For  three  of  thoee.  different  depth  intervale  aqrt  need,  on  one  nation  the  31  meter  depth  interval  lelectad  from  Smith' t 
•ore  included  the  30  meter  interval  of  Moret'i.  on  the  fifth  nation  the  ante  interval  wet  letected 

**  A  multicomponent  model  hae  been  tried  «nth  eeeellent  tuccoM  uam*  cortnuuent  tpectra  for  abaoretion  from  Morel  and  Fneur  <  19TII  Uunt  thoee  functiem  ut*  velvet  of  ft  tat  at  S 
property  choeen  •tvetm|thi.  ,  and e  rjnpe  of  Mapoe  of  the  mooeurad  XIA I  ipectra  could  be  reproduced  and  relative  amount!  of  the  concentration!  of  cMorophylt-lthc  pvnenti.  yetknr 
luonance  and  'other  matenaP  could  be  determined 
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It  was  subdivided  into  5  wavelength  intervals;  (1)  350  to  445 nm,  (2)  445  to  520nm.  (3)  520  to  590nm,  (4)  590  to  660ntn 
and  (5)  660  to  700nm.  and  best  least  square  fits  obtained  to  parabolic  functions.  The  values  obtained  from  the  model  are 
listed  in  Table  2  in  the  “modef  column.  The  model  reproduces  the  values  derived  from  Morel’s  work  with  a  standard 
deviation  of  0.013  as  determined  from  the  ratios  of  the  Model  to  the  Morel  values  (4th  columns  in  Table  2)  and  the  stan¬ 
dard  deviation  of  the  differences  (columns  5)  was  0.0022.  Figures  l(b-f)  show  the  values  plotted  (+>  and  the  analytic  fits 
obtained  (solid  line)  for  each  of  the  5  spectral  regions. 
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Fig.  I.  Comparison  of  values  of  Kw  as  determined  by  means  of  Eq.  (3)  (points  marked  +  )  and  the  fitted  functions  or  'Moder  (shown  i 
the  continuous  line)  Panels  b-f  show  the  S  subregions  that  were  individually  fitted  to  the  points  using  a  least  squares  technique. 


Data  Analysis 

If  one  plots  all  the  available  K  (X)  values  at  each  of  the  68  wavelengths  from  365  to  700  nm  against  the  correspond¬ 
ing  values  of  the  classification  index,  JU490),  one  finds  that  these  scatter  diagrams  provide  strong  evidence  that  there  is  a 
linear  relationship  between  the  two  variables  up  to  some  maximum  value  of  K  (490).  Figure  2.  for  example,  shows  such 
scatter  plots  for  wavelengths  of  400,  440,  520  and  560  nm.  We  have  noticed  that  for  the  shorter  wavelengths  there  is  a 
tendency  for  the  data  to  depart  from  linearity  at  the  larger  values  of  K  (490).  As  a  result  we  have  used  only  data  where 
K  (490)  <  0.160  for  fitting  purposes  when  working  at  wavelengths  below  580  nm.  No  significant  lessening  of  linearity  was 
found  for  longer  wavelengths  when  all  K  (490)  <  0.250  were  included. 

The  linear  equation  is  of  the  form 


K(A)-  /(A)  +  Af(k)  *(490)  (4) 


and  since  ff„(x)  is  a  special  case  of  X(X),  i.e.,  the  lower  limit,  we  see  that  we  can  determine  a  set  of  values  for  AT. (X) 
from  a  knowledge  of  the  slopes.  M( A),  the  intercepts,  /(X),  and  a  value  of  Ku  at  some  reference  wavelength,  say 
K ,  (490)  This  technique  has  been  tried  and  provides  values  very  close  to  those  shown  in  Table  2.  We  have  chosen  to  use 
in  our  analysis,  however,  the  Kw  (X)  values  in  the  columns  labeled  "Moder  (Col.  2)  in  Table  2.  The  slopes  that  were  calcu¬ 
lated  for  each  of  the  68  wavelengths  were  based  on  fitting,  by  least  squares  technique,  a  linear  equation  of  the  form, 

ff(X)-  Af(X)  UU490)-  *.(490)1  +  *„(X)  (5) 
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to  the  data  in  the  scatter  plots.  By  using  the  tabulated  values  or  £„(a)  as  the  intercept  in  the  linear  equation,  we  have  in 
essence  assumed  that  the  £.  ( A)  values  are  correct  The  effect  of  this  was  to  decrease  the  excellence  of  the  lit  slightly  over 
that  which  would  be  obtained  had  we  used  the  data  to  determine  the  values  of  £ „ .  At  most  wavelengths  the  effect  was 
trivial.  Those  few  wavelengths  where  the  effect  seemed  noticeable  were  at  the  extremes  of  the  wavelength  range  where  the 
environmental  and  instrumental  problems  associated  with  field  measurements  were  greater  and  the  number  of  samples  was 
small  or.  in  a  few  instances,  at  wavelengths  where  the  data  was  "noisier'  than  usual. 

Figure  3  is  a  plot  of  the  values  of  the  slopes.  M  (A),  at  each  of  the  68  wavelengths  as  determined  by  the  procedure 
described  above.  A  piecewise  least  squares  fit  was  made  to  these  data  providing  analytic  relationships  that  could  be  used 
for  further  computational  analysis.  Figures  3(b)  to  3(d)  show  the  details  of  the  three  spectral  sections;  (1)  36S-420  nm. 
(2)  420-620nm,  and  (3)  625  to  700nm.  Both  the  data  and  the  mathematical  fit  to  these  data  are  shown.  The  fitting  pro¬ 
cess  not  only  facilitates  the  use  of  the  model  by  the  computer  but  also  provides  additional  smoothing  to  the  overall  data 
set.  essentially  enlarging  its  size  on  a  local  basis  thereby  providing  better  estimates  of  the  slope  than  can  be  obtained  from  a 
regression  at  single  wavelength.  Table  3  lists  the  values  of  the  slopes  as  determined  from  the  individual  fits  to  the  £'(A)  w 
£'(490)  plots  (Col. 3.  headed  "Data")  and  the  values  determined  from  the  fitted  mathematical  functions  (Col. 2.  headed 
"ModeD  The  fourth  column  in  each  panel  is  the  ratio  of  the  value  from  column  2  divided  by  column  3  and  provides  a 
measure  of  the  relative  agreement  between  the  model  and  the  data  used  to  construct  the  model.  The  fifth  column  is  the 
difference  between  the  model  and  the  data  and  is  thus  a  measure  of  the  absolute  agreement.  Slopes  were  calculated  for 
every  5nm.  The  table  lists  values  at  lOnm  intervals. 

Applications 

A  major  application  of  the  tables  is  the  estimation  of  £  (A )  at  a  wavelength  other  than  that  for  which  data  may  be 
available.  For  example,  if  the  value  of  £  (490)  is  known  then  Eq.  (5)  may  be  used  directly  to  compute  £  (A)  Values  for 
Xf (A)  and  £. (A)  are  provided  in  Table  4  at  every  5nm  from  350  to  700nm.  (The  values  of  Af(A)  for  350.  355,  and 
360 nm  contained  in  the  table  were  estimated  by  extrapolation  using  the  fitted  mathematical  function  and  should  be  used 
with  caution.) 

Similarly,  if  the  known  £  is  at  any  wavelength.  A|,  then  the  unknown  £  at  wavelength  Aj  is. 


£(A,)-  [  £  (A i )  —  £„  (A , )  |  +  £.(Aj) 


(6) 
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Table  4.  Values  of  the  slope  function  V/  U  ) 
and  Kw  <A )  for  use  in  Eqs  5  and  6 
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Values  for  S4 (X)  and  KW(K)  should  be  interpolated  when  precise  wavelengths  are  defined.  Equation  (5)  or  (6)  and  Table  4 
are  all  that  are  needed  to  estimate  entire  K  (A )  spectra  if  the  known  value  of  K  indicates  the  water  to  be  oceanic  or 
moderately  clear  coastal  water,  for  example  Jerlov  Type  I  through  III  oceanic  or  Type  1  coastal,  or  if  K (490)  <  0.16  m-1 

Figure  4(a)  is  a  plot  of  the  model  K  (X)  versus  wavelength  with  K  (490)  as  a  parameter.  In  Fig.  4(b)  AT '(A),  i.e. 
K(k)  -  K.  (A),  is  plotted  for  the  same  K (490)  values.  In  Table  5  values  of  K (A)  are  listed  for  each  lOnm  from  350  to 
700nm,  inclusive,  for  a  variety  of  values  of  the  water  classification  index.  K  (490).  from  0.03  to  0.18m-1  Figure  4  and 
Table  5  comprise  a  summary  of  the  results  of  the  model. 


(a) 


(b) 


Figure  4.  Plots  of  the  water  model  with  values  of  KI490I  as  a  parameter  Bottom  curve  K(490>  -  0  022  is  K»U ) 


Figure  5  shows  plots  of  the  model  superimposed  on  IK  (A)  values  calculated  from  the  measured  downwelling  irradi- 
ance  spectra  at  5  of  the  stations  used  in  the  original  data  base.  The  K  (490)  values  range  from  0.033  for  Ftg.  5(a)  to  0  1 16 
for  Fig.  5(e).  The  comparison  between  the  model  and  the  measured  data  provides  examples  of  the  manner  in  which  the 
model  represents  the  data  from  which  it  was  formed.  Figure  6  provides  another  measure  of  the  degree  of  the  fit  between 
the  model  and  the  data.  Figure  6(a)  shows  the  standard  deviation  of  the  ratio  of  the  K  values  calculated  from  the  model 
to  those  obtained  from  the  measured  irradiances.  Figure  6(b)  shows  the  standard  deviation  of  the  differences  between  the 
two.  Figure  7  presents  the  number  of  points  used  at  each  of  the  68  wavelengths  between  365  and  700  in  the  original  data 
set. 
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Fig.  6.  Approximate  standard  deviations  of  comparisons  between  model  and  measured 
data  at  each  of  the  68  wavelengths  examined  between  365  and  700 nm  (a)  Standard 
deviation  of  the  ratio  Kc^)/Km  (X)  where  Kc  is  the  value  calculated  from  the  model  and 
Km  is  ihe  value  determined  from  measured  Ed  ( x  .r )  (b)  Standard  deviation  of  the  dif¬ 
ferences  Kc  IX )  -  Km  (X ) 
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Fig.  7.  Number  of  Km  (X )  samples  used  at  each  wavelength  to  form  the  data  base  for  the  model 
Comparison  with  Jcrlov  water  types 

Ji  is  interesting  to  compare  our  results  with  those  presented  by  Jerlov  (1976).  The  Jerlov  water  types  were  designed 
to  be  used  as  a  general  means  of  classification  the  upper,  usually  well  mixed,  10  meters  of  the  ocean.  The  oceanic  water 
types  I,  I  A,  IB,  II,  and  III  have  a  minimum  K(K)  at  about  475  nm  as  may  be  seen  from  Jerlov’s  Table  XXVII  and  Fig.  70. 
Using  475  nm  as  a  reference  wavelength  we  have  plotted  AC  (X)  vs  K  (475)  using  our  model  and  the  values  of  K  (475)  from 
Jerlov’s  Table  XXVII.  Figure  8  is  a  plot  of  K (X)  from  350  to  600nm  showing  the  manner  in  which  the  attenuation 
changes  for  the  5  Jerlov  oceanic  water  types.  We  observe  that  our  model  shows  the  minimum  attenuation  shifting  to 
longer  wavelengths  as  the  water  turbidity  increases  which  is  consistent  with  most  observations.  We  also  call  attention  to 
the  fact  that  some  of  the  attenuation  values  for  Type  I  oceanic  water  as  given  by  Jerlov  are  less  than  the  values  we  suggest 
for  K.  (X)  In  fact,  at  some  wavelengths  the  values  given  by  Jerlov  for  K  are  less  than  the  values  of  absorption  alone  pub¬ 
lished  by  Morel  and  Prieur  (1977).  We  recommend,  therefore,  that  the  values  of  K  (X)  as  published  by  Jerlov  be  replaced 
by  those  in  Table  6.  The  values  for  Type  I  water  given  in  the  table  are  those  of  pure  sea  water,  Kw .  The  other  tabulated 
values  were  computed  using  our  model  and  the  values  for  K  (475)  given  by  Jerlov  in  his  Table  XXVII.  The  form  of  Table 
6  allows  direct  comparison  with  Jerlov’s  table.  Figure  9  presents  a  comparison  between  the  K  (X)  values  at  the  wavelength 
in  Jerlov's  table  and  our  model  for  oceanic  water  types  IA,  II  and  III.  The  disagreement  between  the  two  is  generally 
greater  in  the  blue.  This  would  be  consistent  with  our  finding  less  evidence  of  yellow  substance  in  the  data  which  we  used 
as  a  basts  for  our  model  than  Jerlov  had  in  the  data  upon  which  he  based  his  work. 
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Fig.  8.  Spectral  diffuse  attenuation  coefficients  for  Jerlov  water  types  Mil. 
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Table  6.  Downward  irradiance  aitenuatiPn  coefficients  K  U)  lor  Jerlov  water  type 
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Fig.  9.  Comparison  between  Jerlov  oceanic  Types  IA,  II  and  III  with  the  present  model 
when  K  (475)  set  equal  to  Jerlov's  tabulated  values. 


Figure  10  is  a  plot  of  K (A)  versus  K( 475)  for  each  25nm  from  350  to  700nm.  This  plot  is  directly  comparable  to 
ure  70  in  Jerlov  (1976).  The  slopes  associated  with  each  of  the  lines  in  the  plot  were  taken  from  the  present  model  and 
listed  in  the  figure. 
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Summary 

The  spectral  characteristics  of  the  diffuse  attenuation  coefficient  of  oceanic  water  behave  in  an  orderly  .  predictable, 
fashion.  The  complete  spectra  may  be  predicted  knowing  the  value  of  the  diffuse  attenuation  coefficient  at  a  single 
wavelength  through  the  use  of  a  model  based  on  the  analysis  of  many  such  spectra.  Regressions  of  the  attenuation 
coefficient  at  any  wavelength  from  365  to  700  nm  against  the  corresponding  value  of  the  coefficient  at  a  second  reference 
wavelength  show  a  linear  relationship  between  the  two  variables  (Eq.  (4),  Fig.  2(a-d))  From  the  slopes  and  intercepts  of 
linear  equations  fitted  to  these  regressions  the  value  of  the  attenuation  of  pure  sea  water.  Kw  ( \ )  may  be  inferred.  Such 
values  were  found  to  be  very  close  to  the  values  obtained  using  Eq.  (3)  with  the  spectral  absorption  and  scattering  values 
published  by  Morel  and  Prieur  (1977)  and  Morel  (1974)  Using  these  latter  values  for  K.  (X)  and  Eq.  (5),  new  slopes. 
.W(x).  were  calculated  which  with  Kw  (A)  form  the  basis  for  the  model  of  spectral  diffuse  attenuation  coefficient.  Table  5 
lists  the  values  of  and  M  at  5nm  intervals  from  350  to  700nm.  These  values  may  be  used  in  Eqs.  (5)  or  (6)  to  pro¬ 
vide  estimates  of  K  (X)  at  any  wavelength,  given  the  value  of  K  at  a  second,  reference  wavelength.  Table  5  lists  values  of 
K  (X)  calculated  as  described  using  selected  values  of  K  (490)  from  0.03  to  0.18m"1  as  a  water  classification  index. 

Jerlov  water  types  may  be  redefined  using  the  model  presented  here.  Table  6  presents  values  of  the  diffuse  attenua¬ 
tion  coefficient  at  the  same  wavelengths  as  published  by  Jeriov  for  oceanic  Types  I.  IA,  IB,  II.  Ill  and  coastal  Type  1.  The 
tabulated  values  at  475  nm  are  the  same  as  those  used  by  Jerlov. 

An  alternative  to  the  use  of  Jerlov  water  types  would  be  the  use  of  a  classification  index  such  as  the  value  of  K  at  a 
selected  reference  wavelength  (e.g.  our  K  (490)).  This  has  the  advantage  that  the  value  of  the  index  provides  a  quantita¬ 
tive  indication  of  the  water  clarity  (i.e.  attenuation  coefficient)  and  allows  the  user  to  classify  the  water  by  selecting  any 
range  of  values  of  the  classification  index  that  may  be  appropriate  to  the  particular  application. 

We  observe  that  the  spectral  nature  of  the  attenuation  coefficient  did  not  appear  to  be  dependent  on  the  geographical 
location  of  the  water  investigated  nor  on  the  depth  of  the  observation  in  the  water  column  providing  the  classification 
index,  K  (490)  was  less  than  approximately  0.16m_l.  We  believe,  therefore,  that  the  model  can  be  considered  as  generally 
applicable  for  water  at  the  surface  or  at  depth  as -long  as  K  (490)  <  0.16m_l. 

References 


Jerlov,  N.G.  (1976).  "Marine  Optics,"  in  Elsevier  Oceanography  Series,  vol.  14,  Elsevier  Scientific  Publishing  Co.,  Amster¬ 
dam,  Oxford,  New  York. 

Morel,  A.  and  L.  Caloumenos  (1971),  "Mesures  D’endairements  Sous  Marins  Flux  De  Photons  Et  Analyse  Spectrale,” 
Centre  De  R echerches  Oceanographiques ,  Report  No.  11.  La  Darse  -  06230  Vi)Jefranche-sur-Mer.  De  Villegranche- 
Sur-Mer,  France. 

Morel.  A.  ( 1974),  “Properties  of  Pure  Water  and  Pure  Sea  Water,"  in  Optical  Aspects  of  Oceanography,  eds  N.G.  Jerlov  and 
B  Steemann  Nielsen,  pp  1-24,  Academic  Press,  New  York,  N  Y. 

Morel,  A.  and  L.  Prieur  (1976),  "Irradiation  Joumaliere  En  Surface  Et  Measure  Des  Edairements  Sous  Marins:  Flux  De 
Photons  Et  Analyse  Spectrale,"  Resultats  De  La  Campagne  Cineca  5  -  J.  Charcot  •  Capricorne  7403,  Report  No  10. 
Sene:  Resultats  des  Campagnes  a  la  Mer,  Groupe  MEDIPROD,  Centre  National  Pour  L'Exploitation  Des  Oceans 
(CNEXO) 

Morel,  A.  and  L.  Prieur  (1977),  "Analysis  of  Variations  in  Ocean  Color,"  Limnol.  Oceanogr..  22(4),  pp.  709-722. 

Okami.  N.,  M.  Kishino,  and  S.  Sugihara  (1978),  "Measurements  of  Spectral  Irradiance  in  the  Seas  Around  the  Japanese 
Islands,"  Technical  Report  No.  2,  The  Physical  Oceanography  Laboratory  of  the  Institue  of  Physical  and  Chemical 
Research.  Tokyo,  Japan. 

Smith,  R.C.  and  K.S  Baker  (1978),  "Optical  Classification  of  Natural  Waters,"  Limnol.  and  Oceanogr  ,  23,  pp.  260-267 
Also  ^ssued  as  University  of  California,  San  Diego,  Scripps  Institution  of  Oceanography,  Visibility  Laboratory.  SIO 

Smith,  R.C.  and  K.S.  Baker  (1981),  "Optical  Properties  of  the  Clearest  Natural  Waters  (200-800nm)."  Appl.  Opt.,  20  pp 
177-184. 

Tyler,  J.E.  and  R.C.  Smith  (1970),  in  Measurement  of  Spectral  Irradiance  Underwater,  pp.  xii  -  103.  Gordon  and  Breach. 
New  York,  London,  Paris. 

Tyler,  J.E.  (1973),  "Data  Report  of  the  Discoverer  Expedition,"  SIO  Ref.  73-16.  pp.  1-1000,  University  of  California.  San 
Diego,  Scripps  Institution  of  Oceanography,  Visibility  Laboratory. 

Acknowledgment 

This  work  was  supported  by  the  Office  of  Naval  Research  under  contract  N  000 14-78- C-0n66.  This  support  is 
gratefully  acknowledged. 


1 78  /  SPI£  Vol.  489  Oeetn  Optics  VII 1 1984) 


